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ABSTRACT 

Aperture  averaging  effects  on  the  measurement  of 
scintillation  and  on  the  temporal -frequency  power  spectrum 
have  been  examined  in  the  marine  boundary  surface  layer 

0 

using  a  He-Ne  (6328  A)  laser  beam  propagating  over  various 
path  lengths  and  under  varying  turbulence  conditions.   Two 
significant  results  have  been  obtained:   (1)  For  scintilla- 
tion measurements,  a  vanishingly  small  receiver  can  best 

be  approximated  by  using  an  aperture  of  diameter 

1/2 
D  .   =  (Z/k)    .   (2)  The  aperture  averaging  correction 

factor  for  an  aperture  diameter  D  >  D  .   decreases  as  the 

c  mm 

level  of  turbulence  increases. 

Data  have  also  been  obtained  supporting  the  theoretical 
prediction  that  an  increase  in  aperture  diameter  reduces 
the  high  frequency  power  content  of  the  temporal-frequency 
power  spectrum. 
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I.   INTRODUCTION 

Scintillation  is  the  term  used  to  describe  the 
atmospherically  induced  intensity  fluctuations  of  a  constant 
power  output  laser  source  as  seen  from  a  distant  detection 
device.   The  precise  manner  in  which  a  detector's  aperture 
size  affects  the  measured  value  of  the  variance  of  the 
intensity  fluctuations  as  a  function  of  distance  from  the 
source  and  lasing  wavelength  has  not  been  satisfactorily 
established.   Several  aperture  averaging  experiments  have 
been  conducted  under  various  turbulence  conditions  to  examine 
this  effect  along  with  an  analysis  of  the  frequency  power 
spectrum  effects  of  aperture  size. 

A.   PROJECT  HISTORY 

Since  1971  an  interdisciplinary  group  of  professors 
principally  from  the  Departments  of  Physics  and  Chemistry, 
Meteorology,  and  Mechanical  Engineering  have  been  investi- 
gating under  a  Navy  research  contract  various  aspects  of 
laser  beam  propagation  over  the  marine  boundary  layer.   Data 
has  been  gathered  thus  far  on  (1)  Modulation  Transfer  Function 
(MTF)  at  wavelengths  of  .4880  urn,  .6328  um,  1.06  ym  and  10.6  ym, 
(2)  scintillation  at  .4880  ym,  .6328  ym  and  10.6  ym,  (3)  extinc- 
tion at  .6328  ym,  (4)  beam  wander  and  (5)  beam  spread. 

The  stated  general  objectives  of  the  project  are:   (1)  to 
develop  an  understanding  of  the  factors  affecting  the  propaga- 
tion of  optical  radiation  through  the  atmosphere,  particularly 
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the  marine  boundary  surface  layer,  the  lowest  30  meters 
above  the  ocean  surface;  (2)  to  develop  the  capability  of 
predicting  optical  system  performance  on  the  basis  of  point 
meteorological  measurements. 

Previous  thought  had  been  that  the  sea  surface  might 
present  a  smooth,  isothermal  layer,  giving  rise  to  little 
turbulence  effect  on  a  propagating  light  beam.   However, 
early  experiments  in  the  project  proved  that  this  was  not 
the  case. 

B.   SCINTILLATION 
1.   Description 

The  twinkling  one  sees  when  observing  distant  light 
sources  is  a  physical  realization  of  scintillation.   This 
appears  as  fluctuation  in  the  output  from  a  small  photo- 
detector  illuminated  by  a  distant  source. 

In  the  general  case  of  finite  beam  propagation, 
measured  intensity  fluctuation  can  be  caused  by  any  or  all 
of  the  following:   (1)  source  output  power  fluctuations, 
(2)  beam  wander,  (3)  detector  motion  through  a  plane 
perpendicular  to  the  optical  axis,  (4)  a  temporal  wavefront 
breakup  within  the  beam  caused  by  small-scale  spatial  and 
temporal  variations  of  the  dielectric  constant  (or  index  of 
refraction)  of  the  lower  atmosphere.   Of  these,  the  most 
critical  for  experimentation  or  information  transmission  are 
the  turbulence  effects  of  scintillation,  beam  wander  and 
beam  spread,  since  they  cannot  be  overcome  by  system 
refinement. 
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2.   The  Effects  of  Turbulence  on  a  Propagating  Wavefront: 
Physical  and  Theoretical  Models 

a.   Kolmogorov  Model  of  Turbulence 

Currently  the  only  well-developed  model  of 
atmospheric  turbulence  is  that  of  Kolmogorov.   It  is  basically 
a  velocity  turbulence  model  that  assumes  homogeneous,  iso- 
tropic turbulence.   The  turbulence  is  assumed  to  be  composed 
of  small  turbulent  eddies,  or  turbules ,  each  of  characteris- 
tic temperature  and  density  (hence  index  of  refraction) ,  as 
illustrated  in  Figure  1.   In  this  homogeneous  turbulence 
model,  often  referred  to  as  "frozen-in  turbulence,"  fluctua- 
tions in  the  dielectric  constant  over  a  vanishingly  small 
volume  element  are  due  to  atmospheric  temperature  variations 
produced  by  the  turbules  being  swept  past  the  reference 
point  by  the  mean  wind. 

Let  n(r,t)  be  the  refractive  index  at  point  r 
in  space  at  time  t.   Then 


(1)        n(r,t)  =  ~     n(r,t')df  +  n.  (r,t',x) 


n(r,t)  =  i-  f    n(r,t' 


lv-' 


=  <n(r,t)>T  +  n1(r,t»,x) 


where  n..(r,t',T)  is  the  deviation  from  the  average  value  at 
point  r  and  time  t. 

The  refractive  index  structure  function  gives 
the  mean  square  variation  in  the  refractive  index  between 
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two  points  in  space: 


(2)  Dn(p)  =  <[n1(r1)  -  n^r^]^  , 


where  p  =  |r..  -  r2 1  . 


The  cross-correlation  function  of  refractive 
index  fluctuations  is  defined  as 


(3)        B 


n(p)  =  i  J    n1(r.1,f)n1(r2,tl)dt' 


=  <n1(£1)n1Cr2)> 


Under  the  assumption  that  homogeneous  turbulence  is  present 
and  the  appropriate  statistical  averages  for  B  (p)  and 
D  (p)  are  stationary  in  time  and  locally  constant  in  space, 
one  finds  that 


(4a)  Dn(p)  =  2[Bn(0)  -  Bn(p)] 


(4b)  Bn(p)  =  y[DnW  -  Dn(p)] 


As  p  becomes  smaller,  correlation  increases  while  the 

structure  function  decreases. 

According  to  the  Kolmogorov  theory  of  turbulence, 

for  I      <<  p  <<  L  ,  that  is,  within  the  "Inertial  Subrange"  of 
oo 

scales , 
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(5)  Dn(p)  =  CN2  p2/3  , 


and  the  corresponding  temperature  structure  function  for 
meteorological  measurements  is 


(6)  DT(p)  =  CT2  p2/3 


iQ   and  L  are  the  inner  and  outer  scales  of  turbulence 
in  the  Inertial  Subrange  (Figure  2) ,  and  CN  and  CT  are  the 
refractive  index  structure  constant  and  temperature  struc- 
ture constant,  respectively.   If  £  (or  L)  represents  the 
diameter  of  a  turbulent  eddy,  I      is  the  scale  size  below 
which  viscous  effects  cause  the  conversion  of  turbulent 
energy  into  heat,  and  subsequent  destruction  of  the  eddy 

ensues.   L  ,  the  outer  scale  size,  is  the  minimum  scale 

o '  ' 

at  which  energy  enters  the  turbulent  system  from  wind  shear 
and  heat  flux. 

For  Kolmogorov  turbulence,  an  analytical  form 
(the  von  Karman  form)  of  the  power  spectrum  (the  three- 
dimensional  Fourier  transform  of  the  cross-correlation 
function) ,  which  illustrates  its  approximate  dependence  on 

I     and  L  ,  is 

o      o 

(7)    *n(K)  =  0.033  Cn2(K2  +  LQ~2)       exp[-(K£Q)2]  . 
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For  K  >>  L  ,  this  reduces  to  the  more  common  form  used  by 
Tatarski , 


(8)       *n(K)  =  0.033  CnV11/3  exp[-K2/Km2]  , 


where  K  =  22L  . 

m  a 
o 


b.   Physical  Models 

(1)  Tatarski 's  Geometrical  Optics  Approach.   The 
physical  mechanism  through  which  turbulence  induces  beam 
breakup  is  by  the  focusing  of  sections  of  a  wavefront 
produced  by  the  inhomogeneities  in  the  index  of  refraction 
over  the  propagation  path.   Tatarski  uses  a  plane  wave 
incident  upon  a  series  of  closely  packed  refractive  spheres 
of  various  sizes.   Using  the  paraxial  approximation,  a 
sphere  acts  as  a  lens  with  focal  length  f  =  * —  ,  where  n, 

is  the  index  of  refraction  perturbation  within  the  sphere. 

For  a  =  1mm,  I     of  the  Inertial  Subrange,  and  n,  =  10    , 

f  =  500m.   For  a  =  lm,  the  lowest  L   limit,  f  =  500  Km. 

'  o      ' 

Thus,  for  most  applications  in  a  turbulent  medium,  a  detector 
will  appear  to  be  closer  to  the  "lens"  than  to  the  focal 
plane.   The  fractional  intensity  change  induced  per  sphere 
will  be  small.   Over  a  long  path,  perturbation  effects  will 
be  multiplicative  and  the  resulting  wavefront  at  the 
receiver  will  be  neither  coherent  nor  of  constant  intensity. 
In  addition,  diffraction  effects,  important  at  distances 
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2 
greater  than  the  Rayleigh  range,  ZR  =  ka  ,  limit  the  ability 

of  the  sphere  to  focus  the  incident  wave.   The  amplitude 

fluctuation  induced  by  a  sphere  at  a  distance  Z  from  a 

detector  is  estimated  [Ref.  1]  as 


a  a    7    2  n,  Z 


From  equation  (9)  it  is  seen  that  the  largest  fluctuations 
are  produced  by  the  smallest  turbules .   By  the  geometrical 
optics  approach  [Ref.  1],  within  the  Rayleigh  range, 
amplitude  fluctuations  increase  as  the  cube  of  the  distance, 
are  dominated  by  the  smallest  turbules,  and  are  independent 
of  wavelength: 


(10)  <^2>   Z    CN2  V7/3  Z* 

o 


However,  this  model  does  not  take  into  account  diffraction 
effects,  important  for  smaller  turbules  that  contribute  most 
toward  amplitude  fluctuations. 

(2)  Diffraction  Effects.   Diffraction  effects 

are  taken  into  account  by  assuming  that  spheres  of  radius 

1/2 
a-  ~  (XZ1)    ,  where  Z1  is  the  distance  from  the  center  of 

the  turbule  to  the  receiver,  contribute  very  weak,  diffracted 

waves  with  no  focusing  effects  at  the  receiver.   Substituting 

1/2 
(XZ)     for  I      in  equation  (10)  ,  the  mean  square  amplitude 

fluctuations  become 
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(11,    <(£,'>  s  c  2  ,7/6  zll/6  (    2   z      2  _ 
o 


6A 
Now  defining  j-   =  X-,  where  X,  is  the  first  order  log- 
amplitude  for  1X2!  <<  1,  Tatarski  has  calculated  the  first 
order  log-amplitude  variance  with  diffraction  effects  con- 
sidered for  plane  waves  as 


(12)  <XX2>  =  0.31  k7/6  CN2  Z11/6 


and  for  spherical  waves  as 


(13)  <XX2>  =  0.13  k7/6  CN2  Z11/6 


2 
both  restricted  to  Z  >>  k£   .   Lutomirski  and  Yura  [Ref.  1] 

state  that  one  would  expect  the  behavior  of  a  beam  wave  to 
lie  between  that  of  the  plane  and  spherical  waves.   There- 
fore, one  may  use  equations  (12)  and  (13)  as  bounds  for 
weak  scintillation.   For  inhomogeneous  turbulence,  equations 
(12)  and  (13)  can  be  expressed  in  the  form 


(14)  <X12>(P)  =  0.56  k7/6  f     (Z-Z')5/6  CN2(Z')dZ»  ,   and 

o 

Z 

(15)  <x12>(S)  =  0.56  k7/6  f     (Z-Z')5/6  (^-)5/6  CN2(Z')dZ'  , 
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where  Z  is  the  distance  from  source  to  receiver  and  Z1 
from  turbule  to  receiver.   For  plane  wave  propagation,  the 
fluctuation  effects  are  weighted  by  the  factor  (Z-Z1)  '  , 
which  has  its  maximum  value  near  the  source  and  approaches 
zero  at  the  receiver.   The  spherical  wave  weighting  factor, 
[(Z-Z')(w— )]    ,  has  its  maximum  value  at  the  midpoint  of 
the  path  and  approaches  zero  at  both  ends. 

c.   Theoretical  Models,  The  Scalar  Wave  Equation 
Numerous  attempts  have  been  made  to  solve  the 
scalar  wave  equation  to  relate  the  electric  field  E(r) 
of  an  optical  wave  passing  through  a  turbulent  medium  to 
the  refractive  index  fluctuations.   Three  approaches  current 
ly  in  use  for  solving  this  problem  are:   (1)  Geometrical 
Optics,  (2)  The  Born  Approximation,  and  (3)  The  Rytov 
Approximation.   Below  is  a  brief  resume  of  the  approaches, 
their  assumptions,  highlights,  methods  of  approach  and 
restrictions.   All  three  methods  start  with  the  scalar  wave 
equation, 


(16)  V2E  +  k2n2E  =  0  . 


The  derivation  of  this  has  assumed  a  medium  with  no  free 
charges  or  current  density,  with  the  permeability  of  free 
space,  that  correlations  between  the  components  of  the 
field  polarized  in  different  directions  can  be  neglected, 
and  that  X   <<    I    . 
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(1)  Geometrical  Optics  Approach 

(a)  Amplitude  fluctuations  are  attributed 
to  the  focusing  and  defocusing  of  the  rays  by  the  curvature 
of  the  turbulent  eddies  along  the  path. 

(b)  Phase  fluctuations  are  calculated  from 
the  changes  in  velocity  as  a  ray  passes  through  regions 

of  different  refractive  index. 

(c)  Diffraction  effects  are  ignored. 

1  /? 

(d)  The  major  restriction  is  (XZ)  '   <<  I   . 

This  reduces  the  applicability  of  this  model  to  extremely 

1/2 
short  ranges.   Note:   When  (XZ)  '   >_  I    ,    diffraction  effects 

become  important. 

(2)  The  Born  Approximation 

(a)  This  is  a  perturbation  technique  to 
incorporate  diffraction  effects  into  the  scalar  wave  equation, 
n  =  <n>  +  n-,  is  substituted  into  the  scalar  wave  equation, 
which  is  solved  for  the  first  order  solution  for  E.   The 
result  is 

(17)       E1  =  J     2k2n1(r»)Eo(r')G(|r-r'|)dV  , 

v 

where  G(|r-r' |)  is  the  free  space  Green's  function. 

(b)  The  perturbation  technique  used  is  valid 
as  long  as  the  amplitude  fluctuations  remain  small. 

(c)  It  is  assumed  that  the  field  at  the 

receiver  is  the  sum  of  the  original  incident  field  and  the 

field  scattered  one  time  from  a  turbule. 
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(d)  It  requires  that  (XZ)  '      >  %    . 

(e)  Restrictions:   In  the  optical  wavelength 

range,  amplitude  fluctuations  often  become  very  large 

compared  to  the  average  value  of  the  amplitude.   Multiple 

scattering  effects  cannot  be  ignored  (but  are) .   References 

2,  3  and  4  contain  details  on  multiple  scattering  effects. 

(3)  The  Rytov  Approximation  (Method  of  Smooth 
Perturbations)" 

(a)  This  is  a  perturbation  technique  applied 

to  a  transformation  of  the  scalar  wave  equation  using  the 

exponential  notation: 

E(r)  =  e^(r)  ,  iKr)  =  X(r)  +  i<t>(r)  . 

to(r)  >  the  "complex  amplitude,"  contains  both  the  amplitude 
and  phase  terms.   Using  this,  the  scalar  wave  equation  takes 
the  form  of  the  non-linear  Ricatti  Equation, 


(18)  V2to  +  Vto-Vto  +  k2(l  +  <5n)  =  0  . 


Now  applying  the  perturbation  technique  to  =  to   +  iK>  where 
ip,  is  the  perturbation  to  the  average  value  of  to,  one 
obtains  the  solution, 

(19)   to1(r)  =  z-jyj  J    G(r-r»)[Vto1.Vto1  +  k26n]EQ  (F')dV'  . 

v 
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The  "Rytov  Approximation"  is  Vi^Vip-  =  |V^- |   =  0,  and 
results  in 

V    f 

(20)     *10(r)  vg^j  J  G(r-r')5n(r')Eo(r')dV  , 

v 

the  "first  order  complex  phase."  A  higher  order  solution 
for  ip,    can  be  obtained  by  iteration. 

(b)  This  approach  includes  multiple  scattering 
effects.   Caveat- -these  effects  are  included  in  an  inflexible 
way  which  does  not  depend  on  the  turbulence  or  other  corre- 
lated parameters  [Refs.  5  and  6]. 

(c)  The  method  requires  that  the  perturbations 
of  the  log-amplitude,  the  real  part  of  the  complex  phase 

[X-l  =  loge(^— )],  and  of  the  phase,  <j>,,  be  small,  i.e., 

X  |  X-.  |    <<    1    and   A  |  <f>    |    <<    1. 

2 

(d)  Restrictions:      X   <<   I   ,    |  Vip  1Q  |      <<    |k   6n|, 

i-£«   1. 

(e)  This  method  includes  both  the  Geometrical 
Optics  approach  and  Born  Approximation  as  special  cases. 

(f)  Theoretical  evidence  indicates  the 

equation  derived  for  amplitude  fluctuations  is  valid  as 

2 
long  as  <x  >  <  0.8  (Saturation  Condition)  [Refs.  7-10]. 

3.   The  Form  of  the  Probability  Distribution  for 
Amplitude  fluctuations 

The  application  of  the  Central  Limit  Theorem  in 

the  derivation  of  the  log-amplitude  equation  leads  to  a 
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prediction  of  a  log-normal  distribution  of  the  amplitude 
fluctuations,  or  equivalently,  a  normal  distribution  for 
log-amplitude.   Some  authors  [Refs.  11  and  12],  taking  a 
different  approach,  argue  for  a  Rayleigh  distribution.   The 
Rayleigh  distribution  arises  when  the  field  at  the  receiver 
is  the  sum  of  the  electric  field  components  from  randomly 
scattered  wavefront  elements.   This  occurs  in  line-of -sight 
propagation  when  the  effects  of  the  turbulent  medium  are 
far  from  the  receiver.   In  line-of-sight  propagation  where 
the  receiver  is  located  within  the  region  in  which  turbulent 
effects  are  being  introduced,  the  field  at  the  receiver  is 
the  result  of  multiplicative  effects  on  a  single  wavefront. 
Strohbehn  [Ref.  2]  shows  that  ip  (t)  ,  representing  the  ampli- 
tude and  phase  terms  of  a  single  wavefront,  is  composed  of 
the  sum  of  a  large  number  of  independent  components,  and 
from  the  Central  Limit  Theorem,  has  a  normal  or  Gaussian 
distribution. 

For  long  enough  paths,  however,  both  the  multiplica- 
tive effects  on  a  single  wavefront  and  additive  effects 
from  different  wavefronts  may  be  important,  and  Strohbehn 
indicates  that  it  is  difficult  to  tell  which  effect  will  be 
predominant.   Empirical  results  thus  far  support  a  log- 
normal  distribution  [Refs.  14,  15  and  16  and  Figure  9]. 
However,  extremely  long-range  experiments  have  not  as  yet 
been  conducted. 
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C.   THE  SATURATION  OF  SCINTILLATION 
1.   Description 

For  plane  and  spherical  waves  (only)  in  homogeneous 
turbulence  and  for  weak  fluctuations, 

a    2 

2  2 

where  a        is  the  normalized  variance  and  a  T   the  log- 

intensity  variance.   In  terms  of  the  log-intensity  variance, 

equations  (12)  and  (13)  become 


(22a)     azl2   =  1.24  k7/6  Z11/6  CN2  ,  (Plane  Wave) 


(22b)     azl2   =  0.52  k7/6  Z11/6  CN2  ,  (Spherical  Wave) 


Empirical  results  [Refs.  14,  32-34]  have  shown  that  when  the 
optically  measured  Gar,    a qj  >    i-s  plotted  against  Tatarski's 

theoretical  values  of  equation  (22),  cj^t    ,  calculated 
from  the  temperature  structure  constant,  the  measured 
scintillations  reach  their  maximum  value  for  tfoT    of  the 
order  of  1  to  3,  then  decrease  slowly  and  tend  to  level  off 
with  increasing  o".j    .   This  approximately  constant  value 
for  Got     >  1  is  called  the  saturation  effect.   Yura 
[Ref.  17]  states  that  the  log-amplitude  covariance  function 
differs  both  quantitatively  and  qualitatively  in  the 
saturation  region  from  that  under  unsaturated  conditions. 
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2 
The  actual  value  of  a   _   at  which  saturation  begins 

has  not  yet  been  strictly  established.   Strohbehn  [Ref.  2] 

indicates  that  saturation  begins  for  a  T   >  0.8,  while 

de  Wolf  [Ref.  18]  and  Lutomirski  and  Yura  [Ref.  1]  maintain 

2 
it  occurs  for  a,.,   >  1. 

2.   Theory 

De  Wolf  has  defined  the  saturation  region  in  terms 

of  an      as 

kL  2 

(23)  1  <  an2    «  ^2-  , 

or  equivalently, 

kL  2 

(24)  1  <  Km1/2  kZ2  Cn2  «  -^-  . 


De  Wolf  [Ref.  18]  improves  upon  the  Rytov  approximation  by 
introducing  a  directional  correction  factor,  cos  [K«p_(Z)  ]  , 
to  allow  rays  being  refracted  by  a  turbule  to  be  deflected 
by  more  than  one  Fresnel  zone  size  off  the  optical  axis. 
This  directional  factor  eliminates  the  focusing  effect  of 
turbules  close  to  the  detection  point  and  introduces  serious 
interference  effects  and  therefore  amplitude  changes. 

De  Wolf  has  calculated  the  log-amplitude  variance  in 
the  saturation  region  based  on  plane  wave  theory  to  be 


(25)      <Xl2>  -  0.4KKJ/3  Z3  CnV1/6  ,  Km  -  fl  , 
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and  derives  the  result  that  the  irradiance  distribution 
remains  log-normal  in  the  saturation  regime.   However,  to 
verify  this  equation,  Km  must  be  determined  empirically, 
which  is  difficult  to  do. 

Other  authors  have  attempted  to  describe  the 
saturation  effect  by  carrying  higher  order  terms  through  the 
wave  equation.   The  efforts  have  been  extremely  complex 
mathematically  and  have  resulted  in  little  success. 

D.   TEMPORAL- FREQUENCY  SPECTRUM  ANALYSIS 
1.   Description 

Like  any  other  time-varying  signal,  a  fluctuating 
voltage  response  from  the  output  of  a  photodetector  possesses 
an  autocorrelation  function.   The  power  spectrum  of  the 
signal  can  be  found  (1)  manually  by  computing  the  Fourier 
transform  of  the  signal's  autocorrelation  function,  or 
(2)  electronically  by  feeding  the  fluctuating  voltage  signal 
into  a  Spectrum  Analyzer,  which  performs  the  transform.   The 
power  spectrum  shows  the  frequencies  composing  the  signal 
and  their  relative  magnitudes.   Current  theories  (described 
below)  indicate  that  the  maximum  frequency  that  will  be 
present  will  be  directly  proportional  to  the  mean  wind  speed 
normal  to  the  propagation  path  and  inversely  proportional 
to  the  square  root  of  the  wavelength  used  and  distance  of 
propagation.   Theories  differ  significantly  on  the  effects 
of  aperture  diameter  on  the  power  spectrum. 
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2.   Theories 


a.   Lutomirski  and  Yura 

It  is  assumed  that  all  refractive  index 

inhomogeneities  of  the  turbulent  medium  move  convectively 

with  mean  wind  speed  v  everywhere  along  the  propagation 

path.   Intensity  fluctuations  will  be  dominated  by  the 

velocity  component  normal  to  the  propagation  path,  v  .   If 

the  transverse  intensity  correlation  length  is  p.,  then  the 

correlation  period  of  the  intensity  fluctuations  at  a  point 

is  of  the  order  of  t   ~  p„/v  ,  giving  a  characteristic 

o    x,  n 

frequency  of 


(26)  f°  =  2^7  • 

1/2 
Using  Lutomirski  and  Yura's  [Ref.  1]  model  of  p„  -  (Z/k)    , 

the  maximum  frequency  cutoff  for  scintillation  effects  at  a 

point  is  given  by 


(2ttXZ) 


The  normalized  photodetector  signal  at  time  t 
is  given  by 

(?^         S-ft)  -  SW  -  <s(*)>  =  J  IC-l,t0  "  <T>  d2p 
(28)    S  (t) <S(t)> /  <T>      a  -1  * 

and  from  the  frozen-in  turbulence  assumption, 
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JCP2  -  V  T,t)   -  <I>    2 

" TT> d  L 2  • 


The  autocorrelation  function  of  S'(t)  is 


(30)   Rs,(x)  =  <s'(t)s'(t+x)>  =  -i-y  rCI(|P1-P2+vix|)d2P1d2P2 


The  power  spectrum  is  the  Fourier  transform  of  R  , (x)  and 


is  given  by 


2ji^V^ 


(31)   *   (a,)  =  -§V-  f    CA^rJl^Fn UJP     ]2  dx 

<I>  u  I  v_  v         ^\/T^ 


n  ^A/l+x' 


Aperture  averaging  effects  on  the  spectrum  are  seen  through 

2v 
the  ooD  term,  ojd  =  -«—  ,  where  D  is  the  aperture  diameter. 

J,  is  a  first-order  Bessel  function  and  C-r  (K)  is  the  two- 
dimensional  power  spectrum  of  the  intensity  fluctuations. 
The  square  bracketed  term  falls  off  rapidly  for  go  >  co^, 
displaying  the  filtering  action  of  the  aperture  diameter. 

When  plotting  a  normalized  power  spectrum, 
$  (f) ,  versus  log  f,  a  change  in  A,  Z,  v  or  D  effects  a 
shift  of  $n(f)  along  the  log  f  axis, 
b.   Clifford  [Ref.  19] 

A  theoretical  solution  for  the  temporal- 
frequency  spectrum  of  scintillation  at  a  vanishingly  small 
aperture  is  given  as  shown  for  plane  and  spherical  waves 
in  Figure  3.   The  corresponding  asymptotic  forms  are  given 
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for  plane  waves   as 


(32)  •(£)    =    0.44    (^£A2)(^-)(1    +    0.27    ft4/3)     ,    ft   «   1    , 

o 

(33)  *(£)    =   1.44    (aM2)(|-)    <T8/3    ,    ft  »   1    , 

o 


and   for   spherical  waves   as 


(34)  <D(f)    =    0.245    (a£A2)(^-)[l    +    0-119    fi4/3]     ,    ft   «   1    , 

o 

(35)  *(£)    =    2.81    (a£A2)(^-)    fi"8/3    ,    ft  »   1    , 

o 

where   ft  =  |-  and   f     =   v    (2ttAZ)"1/2. 
r^  o  n^ 
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II.   NATURE  OF  THE  PROBLEM 

A.   APERTURE  AVERAGING  EFFECTS  ON  SCINTILLATION 
1.   Description 

A  photodetector  placed  at  the  image  point  of  a 
lens  illuminated  with  a  plane  wave  responds  to  the  square 
of  the  integrated  amplitude  over  the  convergent  wavefronts. 
Having  passed  through  a  turbulent  medium,  the  wavefront  of 
incoming  light  for  an  ideal  spherical  or  plane  wave  is  no 
longer  coherent  nor  of  equal  amplitude  across  the  surface 
at  the  detector  aperture.   The  random-sized  turbules 
of  different  refractive  indices  along  the  propagation  path 
have  perturbed  the  wavefront  creating  a  buildup  of  amplitude 
in  some  areas  and  a  cancellation  in  others.   A  transverse 
correlation  distance,  p.,  is  defined  as  the  distance  in  which 
intensity  fluctuations  remain  correlated  on  the  traveling 
wavefront  at  the  point  of  observation,  and  in  general  is 
a  function  of  path  length  Z  and  wavelength  of  the  source. 
If  a  detector  of  diameter  D  is  used  to  measure  the  intensity 
fluctuations  at  a  point  distant  from  the  source,  then  the 
variance  of  the  log-amplitude  measured  will  depend  on  the 
value  of  D  as  compared  to  p..   For  D  >  p»,  the  wavefront 
incident  upon  the  detector  will  not  be  strongly  correlated, 
cancellations  from  the  amplitudes  of  different  parts  of  the 
wavefront  being  out  of  phase  will  occur,  and  intensity 
fluctuations  measured  will  be  less  than  for  a  correlated 
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wavefront,  i.e.,  D  <_  p..   This  observed  decrease  in  the 
variance  of  the  log-intensity  fluctuations  due  to  increasing 
aperture  size  is  called  the  aperture  averaging  effect. 
2.   Theories 

a.   Lutomirski  and  Yura  [Ref.  1] 

For  small  scintillations  and  homogeneous 
turbulence,  the  transverse  correlation  distance  of  the 
intensity  fluctuations  is  given  as 


(36)  p^  s  (Z/k)1/2  =  .399  (AZ)1/2  . 


For  strong  scintillation,  no  reliable  estimate  of  p.  exists 
If  the  aperture  diameter  D  <  p»,  then  the  intensities 
measured  at  all  points  r,  within  the  aperture,  I(r,,t),  are 
correlated  at  time  t,  and  the  photodetector  current  s (t) 
is  proportional  to 


(37)  s(t)  =  ^-  1(0, t)  , 


i.e.,  an  effective  point  detector.   As  D  becomes  greater 
than  po,  s(t)  fluctuates  less  and  for  D  >>  p.,  the  average 
incident  intensity  becomes  constant  and  we  have 


(38)  <s>  =  E2-  <I>  . 


The  signal  variance  is  given  by 


30 


(39)  osZ    =     /    CI(|p1-p2|)d2p1d2p2 


where  C-r(p)  is  the  covariance  of  irradiance.   Lutomirski 

and  Yura  state  that  aperture  averaging  "...is  the  ratio  of 

2 
a        in  (equation  39)  to  what  the  variance  would  have  been 

if  the  intensity  over  the  whole  aperture  varied  as  it  did 

at  a  point."  At  a  point, 


2    ,ttD2 


(40)  as   =  C^-0   Cj(0)  . 


The  aperture  averaging  factor,  A,  is  defined  as: 


Ct  C|  Pi -P?|)   2    2 

(41)      A  '=  f~?  -  tt  I     c;W   *\*\ 


Reference  1  shows  that  if: 

p-  >>  D,  then  A  =  1,  i.e.,  there  is  no  aperture  averaging 
factor; 

Po  2 
p£  «  D,  then  A=  4  (^)  ; 

p«  ~  D,  then  the  value  of  A  requires  a  detailed  knowledge 
of  C-r(p).   Cj(p)  varies  with  source  diameter,  propagation 
slant  path  and  meteorological  conditions.   For  homogeneous 
turbulence  and  plane  wave  propagation,  C-r(p)  is  not  known 
for  strong  scintillation.   It  is  recommended  that 
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(42)  A  =  i _ 

1  ♦  CwELj 


be  used  for  homogeneous  turbulence  when  the  scintillation  is 
small. 

b.   Fried  [Ref.  21] 

Using  the  statistics  of  propagation  of  an 
infinite  plane  wave  over  a  horizontal  path,  Fried  derives 
equation  (37),  then  develops  an  expression  for  the  signal 
variance , 


(43) 


u 

os2    =  2tt  J       pKQ(p,D)  Cj(p)dp  , 


where  K  (p,D)  is  the  area  of  overlap  of  two  circles  of 
diameter  D  whose  centers  are  displaced  a  distance  p,  and 
is  given  by 

Si  {cos"1*)    -    Cfo[l    -    CfoV/2>    ,    P    <   D 
(44)         KQ(p,D)    = 

0,  p    >   D 


Fried  then  derives  the  same  aperture  averaging  factor  as 
given  in  equation  (41)  but  labels  the  factor  0.   He  states 

that  the  normalized  log-amplitude  covariance,  C. (p)/C^  (0) , 

-1/2 
is  a  function  of  p(4Z/k)     only,  for  the  horizontal 

propagation  case,  and  gives  for  an  aperture  diameter 
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normalization  factor, 


(45)  (4Z/k)1/2  =  .798  (AZ)1/2 


Fried's  theoretically  calculated  dependence  of 

1/2 
0  upon  the  normalized  aperture  diameter,  D/(4Z/k)    ,  is 

shown  in  Figure  4  for  various  values  of  C0(0). 

c.   Tatarski  [Ref.  20] 

This  theory  on  aperture  averaging  effects 

gives  the  amplitude  correlation  distance  as  being  "on  the 

1/2 
order  of"  the  first  Fresnel  zone  size,  (AZ)    .   For  this 

model  then,  the  aperture  normalization  factor  is  simply 

1/2 
(XZ)    .   Tatarski  was  the  first  to  define  the  aperture 

„  2 

0  f  R^                 s 
averaging  factor  G(R)  =  )\)r\\    *   wnere  Q(R)  =  t  an<* 

7  ^IUJ  <s>z 

s 
Q(0)  =  — 2_  (  which  Lutomirski  and  Yura  have  labeled  "A" 

<s>z 
and  Fried  "0." 

3.   Published  Experimental  Data 

a.   Homstad  and  Heneghan  [Ref.  22] 

Measurements  were  made  on  a  mesa  near  Boulder, 

Colorado,  over  a  horizontal  (statistically  homogeneous)  path 

at  distances  of  200-1600m  using  a  He-Ne  (6328  A)  laser. 

1/2 
G(R)  vs.  R/(XZ)    was  compared  with  Fried's  theoretical 

curve  of  Figure  4.    G(R)  is  the  aperture  averaging  factor, 

equivalent  to  A  in  equation  (41).   The  result  was  that  good 

agreement  with  theoretical  calculations  was  observed 
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"...in  the  region  of  weak  fluctuations  (i.e.,  at  shorter 

path  lengths) .. .data  in  the  region  of  strong  fluctuations 

1  /? 
...show  a  trend  toward  larger  G(R)  for  R/(XZ)  '      >  1.0,  in 

agreement  with  Kerr"  [Refs.  23  and  24].   This  is  consistent 

with  recent  theoretical  predictions  [Refs.  25  and  26]  and 

with  Kerr's  experimental  work  which  found  the  correlation 

radius  to  be  longer  in  the  region  of  strong  fluctuations 

[Refs.  23,  24,  35  and  36] . 

Homstad  and  Heneghan  point  out  that  a  reason 

that  experimental  data  collected  by  Fried  [Ref.  16]  did  not 

support  his  theoretical  predictions  may  have  been  that  the 

apertures  he  used  for  the  8  km  path  length  were  not  large 

enough  to  extend  beyond  the  correlation  radius  and  hence 

gave  a  relatively  constant  G(R)  as  reported  by  Fried. 

b.  Strohbehn  and  Homstad  [Ref.  27] 

In  an  experiment  conducted  over  land,  good 
agreement  with  Fried' s  model  in  the  region  of  weak  fluctua- 
tions is  confirmed. 

c.  Kerr  [Ref.  23] 

States  that  the  predicted  value  of  the  aperture 

averaging  factor  should  vary  as  D   for  diameters  much 

1/2 
larger  than  (XZ)    ,  where  different  authors  have  assigned 

values  of  2,  7/3  or  3  to  x,  and  that  he  did  not  observe  this 

d.  Previous  NPS  Work  [Ref.  28] 

Schroeder  observed  that  "Significant  aperture 
averaging  was  observed  in  an  enclosed  laboratory  for 
apertures  less  than  the  first  Fresnel  zone  size."   From 
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two  field  experiments  (both  conducted  at  the  same  path 
length  and  with  the  same  wavelength)  Schroeder  stated  that 
the  aperture  averaging  factor  may  be  dependent  on  CN  and 
that  "...a  4.8  mm  aperture  can  approximate  a  vanishingly 
small  receiver  and  should  be  used  whenever  possible."  The 
two  field  experiments  were  conducted  transmitting  across 
Monterey  Bay  with  receiver  and  source  located  on  shore. 
Path  length  was  4.05  km,  98  percent  of  which  was  over  the 
ocean  surface  at  4.9  meters.   Aperture  sizes  used  varied 
from  .3  cm  to  2.3  cm.   By  graphing  "<ao/\(D)  "  cy  .(23mm)  > 
vs.  Aperture  Diameter  D,"  Schroeder's  data  indicated  that 
the  largest  average  difference  in  a.,  values  was  seen  for 
a  4.8  mm  aperture. 

B.   APERTURE  AVERAGING  EFFECTS  ON  THE  TEMPORAL -FREQUENCY 
POWER  SPECTRUM 

1.   Tatarski's  Theoretical  Model 

With  increase  in  aperture  size,  Tatarski  predicts 

a  loss  in  the  high  frequency  components  of  the  normalized 

spectrum.   By  the  addition  of  the  factor 


tVK2  +  47r2f2/v  2 


2J,  (RVK"  +  4ir*fVv 
1 v  n 

RVk2  +  4TTZfZ/v  L 
n 


in  the  kernal  of  the  integral  for  the  plane  wave  model  of 
W(f)  ,  the  spectral  density  of  intensity  fluctuations,  a 
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partial  suppression  of  the  normalized  intensity  fluctuations 

as  seen  by  an  aperture  of  radius  R  is  accounted  for.   For 

v 
f  >>  Jjrrr  ,  this  correction  factor  is  a  very  small  number, 

which  corresponds  to  a  suppression  of  high  frequencies. 
As  the  aperture  size  is  increased,  the  frequency  at  which 
the  suppression  occurs  decreases.   The  low  frequency  portion 
of  the  spectrum  will  remain  unaffected  [Ref.  20]. 

Using  the  frozen-in  turbulence  model,  the  relation- 
ship between  large  spatial  frequencies,  K,  and  large  temporal 
frequencies  is  given  by  to  =  K»v   .   This  implies  that  a 

receiver  of  diameter  D  will  pass  only  temporal  frequencies 
v 

01  10  <  g-  . 

2.   Published  Data 

a.  Hohn  [Ref.  29] 

A  He-Ne  (6328A)  source  was  used  over  path  lengths 
of  4.5  km  and  14.5  km.   Various  aperture  sizes  from  .5  cm 
to  8  cm  were  used.   Wind  speed  values  were  obtained  from 
a  meteorological  center  near  the  test  site.   Normalized 
spectrum  curves  were  deduced  from  the  data  taken  and  the 
results  were  that  "No  systematic  correlations  of  U(f) . . . 
or  f-, /7  were  found  with  changing  aperture  diameter  D  of  the 

photometer."  A  second  result  was  that  the  theoretical 

-1/2 
relation  f  «  Z     did  not  appear  to  hold  empirically. 

b.  Schroeder's  NPS  Work  [Ref.  28] 

An  aperture  averaging  frequency  analysis  sequence 
was  conducted  during  each  of  the  two  field  experiments 
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previously  mentioned.   Schroeder  noted  that  changes  in  wind 
velocity  could  produce  a  different  spectrum  if  two  data 
samples  at  different  aperture  sizes  were  not  conducted 
simultaneously.   General  wind  speed  conditions  were  measured 
at  the  transmitting  location  and  assumed  to  be  the  same 
over  the  propagation  path.   The  spectrum  was  obtained  by 
passing  a  recorded  scintillation  signal  through  a  wave 
analyzer  with  a  bandwidth  of  10  Hz  at  10  Hz  intervals  from 
20  to  100  Hz,  then  at  100  Hz  intervals  from  100  to  1000  Hz. 
The  power  was  determined  by  integrating  the  output  curves 
with  a  planimeter.   Results  were  that  for  f  <  90  Hz,  the 
aperture  averaging  effect  on  the  relative  power  spectrum 
was  "unpredictable,"  while  for  f  >  90  Hz,  the  smaller 
apertures  show  a  larger  contribution  to  relative  power, 
as  predicted  by  Tatarski's  theoretical  model.   A  second 
result  is  that  for  the  same  aperture  size,  the  power  spectrum 
relative  to  the  power  at  20  Hz  measured  on  a  day  when  wind 
speed  was  estimated  to  be  16  knots  over  the  path,  was  overall 
higher  in  magnitude  at  all  frequencies  above  20  Hz  than  for 
the  day  with  wind  speed  estimated  at  one  knot. 

C.   CN0  vs.  CNT 

As  stated  previously,  one  of  the  goals  of  the  project 
is  to  be  able  to  predict,  through  meteorological  measure- 
ments, optical  propagation  conditions.   A  common  indicator 
of  this  is  the  refractive  index  structure  constant,  C^  , 
which  is  related  to  the  log-amplitude  variance  through 
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equations  (12)  and  (13).   CN  can  be  computed  from  the 
temperature  structure  constant,  CT  [equation  (6)],  through 
the  relationship 


6  P 


(46)  CT*  =  (79  x  10  v  iy  )    CNT 


where  CXTT  is  the  refractive  index  structure  constant  pre- 

NT  r 

dieted  from  temperature  fluctuations.   Agreement  between 
the  structure  constant  CNQ  obtained  from  long-path  measure 
ments  of  the  log-intensity  variance  with  CNT  may  be  taken 
as  an  indication  that  the  optical  fluctuations  are  not 
saturated  and  that  the  point-measured  temperature  fluctua- 
tions are  characteristic  of  the  path. 
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III.   EXPERIMENTAL  PROCEDURE 

A.   GENERAL  EQUIPMENT  DESCRIPTION 
1.   Detectors 

To  compensate  for  the  time  dependence  of  the 
variance,  in  measurement  of  effects  of  aperture  size,  a 

United  Detector  Technology  PIN-5DP  photodetector  (Noise 

-13 
Equivalent  Power  5  x  10    watts)  with  fixed  aperture  size 

was  used  as  a  reference.   This  was  operated  without  bias 

voltage  with  the  circuit  shown  in  Figure  5. 

For  the  variable  aperture  detector,  a  TIXL69 
silicon  avalanche  photodiode  (operated  at  its  optimum  S/N 
reverse  bias  voltage  of  159V)  was  used  to  provide  the  high 
sensitivity  at  low  S/N  needed  for  measurement  at  very  small 
aperture  sizes.   The  characteristics  of  the  TIXL69  are 
given  in  Appendix  B. 

For  the  first  series  of  experiments  conducted  in  the 
enclosed  corridor  laboratory,  both  detectors  operated  without 
field  lens  systems.   The  field  lens  images  the  entrance 
aperture  on  the  detector  so  that  the  light  entering  the 
aperture  is  spread  over  the  surface  of  the  detector.   This 
avoids  fluctuations  due  to  the  image  wandering  over  the 
detector  active  area  which  may  vary  in  sensitivity  [Ref.  28]. 
For  the  field  experiments,  the  following  field  lens  systems 
were  installed: 
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(a)  Avalanche  Detector  (Figure  6) .   A  50  cm  focal 
length  objective  lens  with  a  diameter  of  4.9  cm  focused  the 
beam  onto  a  1.5  cm  focal  length,  9.5  mm  microscope  objective 
lens.   Between  the  lenses  a  clear  glass  beam  splitter 
reflected  a  portion  of  the  image  into  an  eyepiece  for  align- 
ment purposes.   The  field  lens  magnification  was  designed 

so  that  the  image  of  the  aperture  filled  the  detector  active 
area  (1.8  x  10"   cm  ) .   A  4.9  cm,  30  A  filter  centered  at 

o 

6328  A  with  a  77  percent  peak  transmission  coefficient 
preceded  the  lens  system,  and  an  18  leaf  iris,  adjustable 
from  1.5  mm  to  42  mm  aperture  diameter,  was  placed  in  front 
of  the  filter. 

(b)  PIN-5DP  Photodiode  Detector  (Figure  7).   This 
field  lens  system  is  fundamentally  the  same  as  that  for  the 
avalanche  detector  with  the  exception  that  a  33.3  cm  focal 
length,  2.54  cm  diameter  objective  lens  and  1.8  cm  focal 

length,  9  mm  diameter  field  lens  were  used  to  cover  the 

2 
detector  active  area  of  .05  cm  . 

2.   Fan-Scan  Mode  of  Laser  Operation 

A  CW  Radiation  5  mWatt  6328  A  laser  operating  in  the 
TEM   mode  was  used  for  the  aperture  averaging  experiments. 
The  NPS  "Fan-Scan"  transmitter  optics  developed  for  ship- 
to-shore  measurements  was  used.   This  system  broadens  the 
circular  beam  cross-section  horizontally  by  reflection  from 
a  variable-radius  cylindrical  mirror,  then  scans  this  "fan" 
vertically  by  reflection  from  an  oscillating  output  mirror. 
Thus  the  narrow  uniform  beam  scans  vertically  twice  each 
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cycle  over  a  detector  within  a  square  field  of  view  of  order 
20  meters  side  at  1500  meters  (variable).   Ship  motion  and 
background  compensation  are  provided  by  measuring  the  dif- 
ference between  maximum  and  minimum  irradiance  values. 
3.   Electronic  Components 

Figure  8  is  a  block  diagram  of  the  electronic 
components  used  in  all  aperture  averaging  and  frequency 
analysis  experiments.   The  voltage  signal  from  each  detector 
package  was  sent  to  a  PAR  113  (Princeton  Applied  Research) 
low  noise  differential  pre-amplif ier  with  maximum  gain  of 
10  K.   The  outputs  from  the  two  PAR's  were  connected  to  a 
single  switch  box,  from  which  the  selected  output  passed 
to  the  demodulator.   This  switch  box  enabled  the  operator 
to  rapidly  record  alternate  scintillation  signals  from  the 
two  detectors.   The  demodulator  stretched  the  voltage  peak 
and  sampled  the  stretched  value,  sampled  the  noise  voltage 
100  microseconds  later  in  between  peaks,  and  sent  the  dif- 
ference to  the  Log  Converter.   The  log-converted  difference 
was  then  sorted  by  the  PIP  into  the  appropriate  voltage 
channel  as  one  count.   The  demodulator  provides  a  trigger 
pulse  for  a  signal  monitoring  oscilloscope  and  the  PIP. 

The  calibration  box  provides  calibration  voltages 
of  .5,  1,  2,  4  and  8  volts  to  the  Log  Voltmeter/ Converter , 
and  from  there  into  the  appropriate  voltage  channel  on  the 
PIP  400  Multichannel  Pulse  Height  Analyzer  (MPHA) .   A 
readout  of  the  counts  in  each  channel  during  calibration 
is  used  by  the  HP9810A  Calculator  and  Plotter  to  scale  the 
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voltage.   The  PIP  400  MPHA  records  the  number  of  times  a 
particular  voltage  value  (i.e.,  log-intensity)  is  measured 
from  the  fluctuating  signal  during  a  run.   This  information 
is  displayed  durirtg  collection  and  can  be  stored  indefinitely 
for  subsequent  printout  and  analysis  as  long  as  no  other 
data  is  recorded  further  on  that  channel.   The  scintillation 
signal  for  frequency  analysis  was  recorded  from  the  output 
of  the  Log  Voltmeter/Converter. 

4.   Shipboard  Detector  Stabilization 

The  avalanche  and  PIN  photodiode  detectors  were 
mounted  on  a  feedback-controlled  gyrostabilized  platform. 
The  system  is  designed  to  compensate  for  ship  motion  and 
keep  the  critical  detector  alignment  stable.   Currently  this 
system's  performance  is  adequate  at  short  ranges  (less  than 
2000  m)  and  fairly  calm  sea  states,  but  is  degraded  rapidly 
by  increased  distances  or  moderate  sea  states.   In  these 
instances  the  detector  (gyroscope)  alignment  must  be  con- 
trolled manually  by  keeping  the  laser  signal  image  in  the 
center  of  two  crosshairs  of  an  eyepiece  mounted  on  the 
gyro  support.   Gyroscope  nutation,  inability  to  keep  the 
detectors  aligned  properly,  and  movements  of  the  ship  all 
contribute  errors  to  a  measured  value  of  signal  variance. 

B.   APERTURE  AVERAGING  EXPERIMENTS 

1.   Propagation  Paths  and  Turbulence  Conditions  Present 
a.   Corridor  Laboratory 

Aperture  averaging  experiments  were  conducted 
in  a  long  basement  corridor  by  locating  a  5  milliwatt  He-Ne 
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(6328  A)  laser  at  one  end  of  the  corridor,  propagating  the 
beam  the  length  of  it,  reflecting  it  back  down  the  length 
of  the  corridor  using  a  plane  mirror  and  measuring  the 
scintillation  near  the  location  of  the  source.   Total  path 
length  was  280  m.   The  avalanche  and  PIN  photodiode 
detectors  were  clamped  to  a  vertical  support.   Once  aligned, 
the  detectors  were  not  subject  to  any  of  the  inherent  motions 
from  platform  instability  which  are  present  when  conducting 
the  same  experiment  on  board  a  ship.   The  experiments  were 
conducted  during  early  morning  hours  (with  the  exception 
of  that  on  7  April)  to  eliminate  rapid  variances  in  turbu- 
lence conditions  due  to  persons  opening  and  closing  exterior 
doors.   Turbulence  was  provided  by  ventilation  system  openings 
and  building  heating  ducts  located  near  the  propagation  path. 
Earlier  measurements  on  the  nature  of  the  turbulence  in 
this  corridor  indicated  that  the  turbulence  present  did  not 
represent  homogeneous  (Kolmogorov)  turbulence  conditions. 
However,  the  non-homogeneous  turbulence  did  produce  an 
intensity  distribution  that  proved  to  be  log-normal  as 
predicted  from  theory.   A  sample  log-intensity  distribution 
taken  during  a  corridor  experiment  is  shown  in  Figure -9. 
b.   Over  Land 

To  complete  an  aperture  averaging  sequence  in 
homogeneous  (Kolmogorov)  turbulence  before  going  to  sea, 
where  the  ship's  motion  can  affect  the  measured  signal 
variance,  an  aperture  averaging  sequence  was  conducted  on 
the  Naval  Postgraduate  School  Annex  grounds  on  which  the 
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Fleet  Numerical  Weather  Central  is  located.   The  site 
selected  to  conduct  the  experiment  consisted  of  a  level 
terrain,  three  fourths  of  which  was  a  blacktop  covered  road, 
the  remainder  a  level  dirt  road.   There  were  buildings  and 
trees  within  50  m  of  the  propagation  path;  however,  no 
forced  draft  blowers,  heating  systems  or  machinery  was 
present  to  disturb  the  apparent  homogeneity  of  the  path. 
Winds  were  negligible;  hence  no  eddy  effects  from  winds 
curling  around  buildings  or  trees  were  present.   The  5  mW 

o 

He-Ne  (6328  A)  source  was  mounted  on  a  large  telescope 
frame  in  a  small,  open-end  trailer.   It  transmitted  a 
Fan-Scan  signal  over  a  path  length  of  336  m  to  the  avalanche 
and  PIN  photodiode  detectors  mounted  on  a  metal  plate 
affixed  to  a  manually  trained  telescope.   The  equipment 
setup  of  Figure  8  was  used  and  a  frequency  analysis  experi- 
ment was  conducted  in  conjunction  with  the  aperture  averaging. 
Again,  intensity  distributions  over  the  path  proved  to  be 
log-normal. 

c.   Over  Water  (Land  to  Land) 

Two  experiments  were  conducted  over  different 
path  lengths  and  under  different  turbulence  conditions. 

Experiment  I:   Coast  Guard  Wharf  to  Spanagel 
Roof.   The  optical  path,  of  length  2335  m,  was  highly 
irregular.   The  propagation  path  was  slanted  downwards, 
going  from  the  roof  at  an  altitude  of  about  75  m  above  sea 
level  down  to  the  Mobile  Optical  Research  Laboratory  at  the 
Coast  Guard  Station  parking  lot  at  about  10  m  above  sea  level. 
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About  one  third  of  the  path  was  over  water,  the  remainder 
over  several  commercial  and  residential  structures.   During 
the  course  of  the  experiment  there  were  intermittent  rain 
squalls.   Data  revealed  that  this  entire  sequence  was 
conducted  under  saturation  conditions  for  a.-r  . 

Experiment  II:   Monterey  Holiday  Inn  to  Hopkins 
Marine  Station.   The  second  experiment,  conducted  over  a 
propagation  path  length  of  4230  m,  was  97  percent  through 
the  marine  boundary  layer.   The  laser  beam  was  transmitted 
over  a  horizontal  path  at  about  10  m  above  sea  level.   The 
sky  was  clear  on  11  March  with  average  wind  speeds  during 
the  experiment  of  about  7  m/sec.   The  He-Ne  (6328  A)  laser 
source  was  located  in  the  rear  of  the  open-end  trailer  and 
was  operated  in  the  Fan-Scan  mode.   The  trailer  was  located 
at  the  Monterey  Holiday  Inn  parking  lot,  about  100  m  from 
the  water's  edge.   The  terrain  between  the  ocean  and  the 
trailer  consisted  of  a  low,  sandy  incline  sloping  downwards 
from  the  trailer's  position  towards  the  ocean.   The  detection 
devices  and  processing  system  were  located  in  the  bus  (Mobile 
Optical  Research  Laboratory)  stationed  on  the  Hopkins  Marine 
Station  grounds  about  10  m  above  the  ocean  surface  and  20  m 
from  the  water's  edge  with  an  unobstructed  view  of  the 
signal  source. 

Experiments  measuring  MTF  being  conducted 
simultaneously  with  aperture  averaging  precluded  the  PIN 
photodiode  detector's  being  available  for  use  as  the  fixed 
aperture  system  on  this  occasion.   For  this  experiment  a 
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HgCdTe  detector  with  an  aperture  size  of  23.6  mm  was  used 
as  the  fixed  aperture  system.   The  HgCdTe  detector  measured 

scintillation  from  a  (10.6  urn)  CO  laser.   It  is  assumed 

2 
that  CN  measured  at  10.6  urn  varied  temporally  in  the  same 

way  as  measured  at  .6328  urn,  as  indicated  by  the  consistency 

2 
of  multiwavelength  measurements  of  C„   . 

d.   Over  Water  (Land  to  Ship) 

The  He-Ne  Fan-Scan  transmitter  on  the  bus  was 
located  at  Pt.  Pinos  (Figure  10)  as  near  to  the  shoreline 
as  possible  (about  50  m  away) ,  with  the  detectors  and  data 
processing  system  on  board  the  Research  Vessel  ACANIA.   Only 
those  runs  were  recorded  in  which  detector  alignment  was 
maintained  for  periods  of  at  least,  two  to  three  minutes. 
Path  lengths  used  during  these  runs  were  1442  m,  1550  m, 
and  1625  m;  the  maximum  range  was  limited  by  the  failure 
of  the  tracking  system  to  lock-on  automatically  at  ranges 
greater  than  2000  m.   Path  lengths  were  determined  by  base- 
line triangulation  methods  on  shore  using  surveyor's 
transits.   The  optical  paths  were  about  98  percent  over  the 
ocean  surface. 

2.   Data  Reduction 

Data  reduction,  i.e.,  determining  the  refractive 
index  structure  constant  and  the  log- intensity  variance 
from  a  fluctuating  voltage  signal,  was  accomplished  within 
minutes  of  the  completion  of  a  data  collection  run  through 
the  use  of  the  Victoreen  PIP  400  Multichannel  Pulse  Height 
Analyzer  and  the  Hewlett  Packard  9810A  Calculator  and 
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Plotter.   As  described  previously,  the  PIP  400  stored  the 
log-intensity  distribution  of  the  irradiance  measured. 
An  interface  was  developed  between  the  PIP  400  and  the 
HP9810A  to  allow  the  calculator/plotter  to  print  out  the 
log- intensity  distribution  and  determine  CN  and  a.,  from 
it.   The  algorithm  the  computer  followed  [Ref.  31]  was  to 
apply  a  best- fit  Gaussian  curve  to  the  distribution  to 
obtain  o"»,,  and  then  determine  CN  through  the  relation 


(47)  CM   = 


k//12    Zli/12    (b476)1/2    ' 


An  example  of  the  printout  is  shown  in  Figure  9.   The  three 
markers  near  channels  60,  116  and  173  are  calibration 
points  of  .5  V,  2  V  and  8  V  respectively,  which  are  equally 
spaced  on  the  log  scale. 

3.   Aperture  Averaging  Factor  G(D) 
a.   Description 

To  determine  the  aperture  size  that  best 
approximates  a  point  receiver,  i.e.,  that  size  at  which 
a II   reacnes  a  maximum  for  a  given  path  length  and  C^ 
value  and  does  not  decrease  significantly  for  further 
reduction  of  aperture  size,  the  following  aperture  averaging 
factor  was  used: 
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urn       rrm   -  ^I"W]       V^V«]       cN[V(D0>t)I 

(48}      G(D)  =  cNi^cu0,t0jj   *  CNlV(l),tJJ        cNmu,tjj     • 


Cvj[V(D  , t  )]  is  the  value  of  CN  measured  using  the  variable 

aperture  detector  at  the  maximum  diameter  D  and  tne 

time  t  .   Similarly,  F  denotes  the  fixed  aperture  detector 

which  is  used  to  measure  the  change  in  CN  due  to  variables 

other  than  a  change  in  aperture  size.   The  first  factor 

CN[V(D0,to)]/CN[F(D0,tQ)],  compares  the  initial  CN  value 

measured  (at  time  t^)  with  the  variable  diameter  set  at  D^ 
v         o  o 

to  that  measured  by  the  fixed  diameter  at  time  t  . 
CN[F(D0,t)]/CN[V(D,t)]  reflects  the  change  in  that  ratio 
produced  by  a  reduction  in  aperture  diameter  from  D  to  D. 
If  the  two  apertures  (fixed  and  variable)  measure  a»j  simul- 
taneously at  times  t  and  t,  and  the  variable  aperture  size 
is  not  changed  between  the  two  measurements,  then  G(D) 
will  ideally  equal  1.00.   This  was  tested  in  conjunction 
with  two  corridor  experiments  conducted  on  10  February 
and  5  March  19  76.   The  comparison  runs  produced  G(D)  values 
of  .9818  and  .9981,  indicating  an  average  error  of  about 
one  percent  for  G(D)  values  at  different  aperture  sizes 
due  to  experimental  variances. 

b.   Reason  for  Use;  Relation  to  Fried's  0(D) 
By  starting  with  the  largest  experimental 
variable  aperture  diameter,  D  ,  for  which  G(D  )  =  1.00, 
then  progressively  reducing  it  until  the  minimum  variable 
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aperture  D  -n  is  reached,  the  aperture  size  that  best  approxi- 
mates a  point  detector  can  be  determined  by  noting  at  which 
aperture  diameter  the  lowest  value  of  G(D)  occurs.   Empirical 
results  (Section  IV. A.)  show  that  G(D)  does  not  remain  linear 

from  D  down  to  the  smallest  possible  aperture  size,  D    , 

o  r  '      mm* 

but  levels  off  and  reaches  a  minimum  value,  [G(D)]  .  ,  well 

before  it. 

Fried  [Ref.  21]  defined  an  aperture  averaging 

factor   0(D)   [equivalent  to  A  in  equation  (41)]  based  upon 

a  vanishingly  small  receiver  (D  ■*•   0  mm).   Fried's  theory 

states  that  all  apertures  of  diameter  D  >  0  will  measure  a 

signal  variance  of  smaller  magnitude  than  would  be  measured 

using  a  vanishingly  small  detector.   When  using  an  aperture 

of  diameter  D,  the  measured  variance  must  be  multiplied  by 

the  correction  factor   grrn   to  allow  for  aperture  averaging. 

TG(D)]  •   identifies  the  aperture  diameter,  D  .  ,  that  best 
L  ^  *  Jmm  r  '   mm* 

approximates  a  vanishingly  small  receiver.   G(D)  is  related 
to  Fried's  aperture  averaging  factor  in  the  following  manner: 


r„Q1         GCD)  CN[V(D0,t)]         CN[V(D,t)]  CN[V(D0,t)]  x 

C      J    ™^In     ^NLV(.1>.*JJ       '    CNlVtJJnin,tjJ-   CnlH\in,t)]-    5W 


C.   APERTURE  AVERAGING  EFFECTS  ON  TEMPORAL -FREQUENCY 
SPECTRUM  ANALYSIS 

1.   Experimental  Conditions 

a.   (7-9)  April  1976,  Basement  Corridor 

The  electronic  equipment  setup  used  in  aperture 

averaging  sequences  is  shown  in  Figure  8  and  the  procedure 
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followed  in  obtaining  and  processing  signals  on  this 
date  is  described  in  Sections  III.B.l.a,  IV. A. 1,  and 
IV. A. 2. a. (3) .   For  frequency  analysis,  the  scintillation 
signal  was  recorded  from  the  output  of  the  Log  Voltmeter/ 
Converter  with  a  Precision  Instruments,  Model  No.  PI -6204, 
four-channel  recorder  at  37.5  ips.   Performance  characteris- 
tics of  the  recorder  are  contained  in  Appendix  C.   Signals 
were  recorded  for  a  period  of  about  75  seconds  each  at 
aperture  sizes  of  40,  30,  20,  15,  10,  7,  5,  3,  2  and  1.5 
millimeters . 

b.  20  April  19  76,  Navy  Annex 

Experimental  conditions  are  as  described  in 
Sections  III.B.l.b,  IV. A. 1,  and  IV.A.2.C.   Frequency  analysis 
data  was  recorded  as  described  above  for  the  (7-9)  April  run. 
Aperture  diameters  examined  were  of  40,  30,  20,  15,  10,  7, 
5  and  3  millimeters.   Power  to  operate  the  equipment  on 
the  bus  was  provided  by  a  portable  gasoline-powered  motor 
generator. 

c.  (29-30)  April  1976,  R/V  ACANIA 
Experimental  conditions  on  this  date  and 

procedures  followed  for  aperture  averaging  are  described 
in  Sections  III.B.l.d,  IV. A. 1,  and  IV.A.3.(d  and  e) .   The 
signal  to  be  frequency  analyzed  was  recorded  from  the 
output  of  the  Log  Voltmeter/Converter  (Figure  8)  on  a 
Sanborn  HP- 3914,  14-channel  Analog  Tape  Recorder,  operated 
by  the  meteorology  team  recording  micrometeorological  data 
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for  CM  comparison.   Operating  characteristics  of  the 
HP-3914  Recorder  are  included  in  Appendix  C.   Time  reference 
signals  were  recorded  on  the  tape  for  identification. 
2.   Data  Reduction 

The  recorded  scintillation  signal  was  frequency 
analyzed  by  playing  back  the  signal  into  an  EMR  Schlumberger 
Instruments  1510  Digital  Spectrum  Analyzer.   The  spectrum 
analyzer  was  operated  in  the  continuous  mode  over  a  frequency 
range  of  10  Hz  to  2.56  KHz.   The  signal  was  sampled  by  the 
spectrum  analyzer  for  a  period  of  25.6  seconds.   The  resulting 
spectrum  was  plotted  on  a  Hewlett  Packard  Moseley  7035 B  X-Y 
Recorder.   Frequency  was  plotted  linearly  along  the  X-axis 
from  42  Hz  to  2.56  KHz.   The  Y-axis  was  a  0  to  -60  dB  linear 
scale. 

For  evaluation  of  the  aperture  averaging  effects  on 
the  spectral  density,  relative  power  spectra  were  obtained 
directly  by  aligning  the  -60  dB  line  of  the  spectrum  of  each 
different  aperture  to  a  common  new  reference  axis  and 
plotting  each  spectrum  manually.   Points  of  local  minima 
were  plotted  on  each  spectrum,  thereby  avoiding  abnormal 
magnitudes  due  to  noise  spikes. 
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IV.   EXPERIMENTAL  RESULTS 

A.   APERTURE  AVERAGING  EFFECTS  ON  SCINTILLATION 

1.  Commonalities 

In  addition  to  the  electronic  equipment  setup 
(Figure  8),  detectors,  and  laser  source  previously  described, 
the  following  items  were  common  to  all  aperture  averaging 
experiments:   (1)  The  Hewlett  Packard  9810A  Calculator  and 
Plotter  was  used  for  immediate  data  reduction.   (2)  Aperture 
sizes  for  the  variable  aperture  system  were  set  by  inserting 
a  micrometer  caliper  set  at  the  proper  diameter  size  within 
the  aperture  and  closing  the  iris  upon  it.   Both  horizontal 
and  vertical  dimensions  were  checked.   The  aperture  size  was 
rechecked  at  the  end  of  the  run. 

2.  Preliminary  Investigations 
a.   Corridor  Experiments 

(1)  2235,  7  February  1976  -  0515,  8  February  1976 
This  initial  aperture  averaging  sequence  was 
conducted  to  examine  the  validity  of  the  assertions  made 
by  Schroeder  [Ref.  28]  that  a  4.8  mm  aperture  diameter  best 
approximates  a  vanishingly  small  receiver.   For  this  experi- 
ment, the  laser  was  operated  in  the  Scan  mode  only  and  a 
neutral  density  filter  of  .5  was  used  to  decrease  the  signal 
strength.   The  field  lens  detector  systems  of  Figures  6  and 
7  had  not  been  constructed  at  this  time.   Each  detector 
consisted  of  the  photodiode  placed  at  the  focal  plane  of  an 
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objective  lens  in  a  narrow  cylinder  of  diameter  22.6  mm  with 

a  10  X  filter  centered  at  6328  A  at  the  cylinder  front.   This 

made  the  largest  variable  aperture  size  available  22.6  mm. 

Aperture  sizes  of  *22.6,  20,  17,  14,  12,  10,  7,  6,  5,  4,  3 

and  2  millimeters  were  used.   G(D)  values  calculated  from 

the  measured  log-intensity  distribution  standard  deviations, 

a« y,  are  given  in  Table  I  with  the  plotted  MG(D)  vs.  Aperture 

Diameter,  DM  shown  in  Figure  11.   The  last  column  of  Table  I 

1         G(D)       .  - 

gives  q(d)       =  rQ(u)  i »  tne  correction  factor  for  aperture 

*■    *■  ^-"min 
averaging. 

In  Figure  11,  the  curve  bottoms  out  at  an 
aperture  size  of  about  5  mm  and  the  corresponding  value  of 
G(D)  does  not  significantly  decrease  from  its  value  at 
that  point  as  the  aperture  diameter  is  further  reduced. 

(2)  1655-1820  and  2100-2240,  10  February  1976 

It  was  noted  from  the  experiment  of 
7  February  that  although  the  log-intensity  distributions 
were  for  the  most  part  Gaussian,  on  several  of  the  larger 
aperture  runs  double  peaks  occurred  near  the  top  of  the 
log-intensity  curve.   It  was  believed  that  this  was  due  to 
the  laser  being  in  the  Scan  mode  only  so  that  the  detector 
would  see  one  average  intensity  peak  as  the  beam  swept  up 
and  a  slightly  different  one  as  the  beam  swept  down.   In  an 
attempt  to  alleviate  this,  a  second  aperture  averaging 
sequence  was  conducted  with  the  laser  in  the  Fan-Scan  mode. 
By  spreading  the  beam  out  and  scanning  it,  it  was  felt  that 
a  more  uniform  intensity  might  be  seen  by  the  detector,  and 
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the  double  peaks  eliminated.  A  Fan-Scan  divergence  setting 
of  seven  milliradians  was  used  and  a  neutral  density  filter 
of  optical  density  .5  was  inserted  in  the  beam  to  reduce  the 

signal  intensity  at  the  detector.   With  the  beam  reflected 

2 
off  a  625  cm  mirror  at  a  distance  of  140  m,  the  estimated 

-9 
signal  strength  at  the  position  of  the  detector  is  1.7  x  10 

2 
Watts/cm  ,  which  is  of  the  order  of  the  signal  intensity 

for  the  at-sea  experiments,  described  later,  at  a  path 
length  of  1600  m. 

Aperture  sizes  used  were  identical  to  those 
used  during  the  7  February  run.   Log  Intensity  distributions 
on  this  day  revealed  that  the  double  peaks  were  eliminated 
in  all  but  one  distribution  by  use  of  the  Fan-Scan,  and 
this  mode  was  therefore  adopted  as  the  standard  laser  opera- 
tion mode  for  future  aperture  averaging  sequences.   Figure 
12  shows  MG(D)  vs.  Aperture  Diameter,  DM  for  this  experiment. 
Although  this  experiment  was  conducted  identically  to  that 
on  7  February  with  the  exception  of  the  Fan-Scan  operation, 
three  significant  differences  stand  out  in  the  results. 
First,  the  curve  bottoms  out  at  the  6  mm  aperture  size  vice 
the  5  mm  size  of  7  February.   Secondly,  G(D)  values  at  all 
aperture  sizes  are  in  general  larger  than  they  were  on 
7  February,  even  though  the  same  value  of  D  (22.6  mm)  was 
used  to  calculate  G(D)  values.   Thirdly,  the  shape  of  the 
curve  following  the  bottoming  out  at  6  mm  is  different  from 
that  seen  on  7  February.   In  comparing  the  data,  the 
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average  a. ,.  [F(D  ,t)  ]  is  seen  to  be  significantly  different 
on  the  two  days,  implying  a  larger  CN.   On  7  February 
<an  [F(DQ,t)]>  =  .330,  while  on  10  February  <o ^ [F(D  ,t) ] >  =.428 
An  increasing  <o'of>  seems  to  increase  G(D)  values,  as  can  be 
seen  by  comparing  Table  II  with  Table  I. 
(3)  (7-9)  April  1976 

An  aperture  averaging  frequency  analysis 
experiment  was  conducted.   The  aperture  averaging  results 
are  included  for  comparison  with  runs  conducted  on  7  February 

o 

and  10  February.   The  same  He-Ne  (6328  A)  laser  was 

operated  in  the  Fan-Scan  mode  with  a  divergence  setting  of 

seven  milliradians  and  a  neutral  density  filter  of  .5  was 

inserted  in  the  beam. 

Prior  to  this  run  the  largest  aperture  size 

available  on  the  variable  aperture  detector  was  22.6  mm. 

To  increase  the  largest  aperture  size  available  and  improve 

signal  detection  capability,  the  field  lens  systems  of 

Figures  6  and  7  had  been  constructed  and  were  in  use  for 

this  experiment.   The  significance  here  is  that  with  a 

larger  D  available,  one  would  expect  to  see  a  smaller 

[G(D)]  •   for  this  sequence.   This  proved  accurate,  as 

shown  in  Figure  13.   A  much  lower  fG (D")  ]  •   resulted  using 

L    Jmin 

the  40  mm  aperture  as  D  vice  22.6  mm.   This  curve  appeared 
to  follow  the  same  basic  pattern  as  those  for  7  February 
and  10  February  with  the  exception  that  this  took  on  a 
different  shape  from  either  of  the  previous  two  for  apertures 
less  than  the  bottoming  out  diameter  of  about  5  mm.   In 
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order  to  see  better  how  this  run  compared  with  the  previous 

two,  G(D)  values  were  recomputed  using  the  20  mm  diameter 

as  D  .   The  result  is  shown  in  Figure  14.   [G CD")  ]     for 

(7-9)  April  is  still  well  below  that  of  7  and  10  February. 

The  average  a--,  for  the  fixed  aperture  system  for  this  run 

was  <cr. -j.  [F(D  ,t)  ]>  =  .08  which  is  much  lower  than  the  average 

values  for  7  and  10  February  as  is  [G(D)1  .  .   Table  III 

L    Jmin 

gives  a  comparison  of  the  average  signal  variance  with 
[G(D)j  •   for  the  three  experiments  conducted  in  the  same 
manner,  with  the  same  equipment,  but  with  different  turbu- 
lence conditions  present  as  indicated  by  <^j,I>  . 


TABLE 

III 

Date 

i 

<0i 

[F(D  , 

t)] 

>     [G(D)1  • 

L  v  '  Jmm 

^liy 

A[G(D)]  . 

L  *■  *  Jmm 

10  Feb 

.428 

.647 

7  Feb 

.330 

.468 

.098 

.179 

(7-9) 

Apr 

.084 

.223 

.246 

.245 

It  appears  that  [G(D)1  •   increases  with  <onT>   for  low  values 

rr  L  v  J  Jmm  ill 

of  <aoT>  ,  and  that  in  general,  for  a  given  D  ,  the  value  of 
[G(D)]  .   is  dependent  upon  the  level  of  turbulence.   The 
higher  the  turbulence  level,  the  higher  the  [G(D)]min  value, 
b.   Aperture  Diameter  Normalization 

Examination  of  Figure  14  shows  that  the  shapes 
and  slopes  of  the  plots  are  similar  from  the  largest  aper- 
ture size  down  to  about  5  mm  where  three  different  forms  of 
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the  curve  then  emerge.   Following  previous  theoretical 

models  [Refs.  1,  20  and  21],  it  is  expected  that  D    ,  the 

r  mm 

aperture  best  approximating  the  "vanishingly  small"  receiver, 
should  be  directly  related  to  the  intensity  transverse 
correlation  length,  p..   This  quantity  is  given  in  the 
three  models  as : 

(1)  (AZ)1/2  ,  [Tatarski,  Ref.  20] 

(2)  (^)1/2   ,  [Fried,  Ref.  21] 

(3)  C^)1/2    ,  [Lutomirski  and  Yura,  Ref.  1]. 

Aperture  diameters  for  the  three  corridor  experiments  were 

scaled  according  to  the  three  models,  the  results  of  which 

are  shown  in  Figures  15,  16  and  17,  respectively.   The 

three  curves  bottom  out  at  approximately  the  value  of 

Z  I/2 
D  =  (V)    predicted  by  Lutomirski  and  Yura. 

c.   Over  Land  Experiment  (2128-2347,  20  April  1976) 

Aperture  diameters  examined  were  40,  30,  20, 

15,  10,  7,  5  and  3  millimeters.   The  Fan-Scan  was  set  at  a 

divergence  of  7  milliradians  and  no  neutral  density  filter 

was  used  to  decrease  the  signal  intensity.   The  laser  was 

situated  about  2  meters  and  the  detectors  2.5  meters  off 

the  ground.   The  sky  was  clear,  temperature  averaged  10°C 

and  winds  were  negligible  (0-2  m/sec) . 

MG(D)  vs.  D  , ,-"  was  plotted  to  locate  the 

(Z/k)i/Z 
bottoming  out  point.   Although  able  to  test  only  a  few 

aperture  sizes,  previous  runs  conducted  outside  the  corridor 
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laboratory  (discussed  in  section  d  below)  indicated  that  the 

1/2 
(Z/k)    normalization  factor  was  accurate  in  predicting 

1/2 
the  bottoming  out  point  of  the  graph  at  D/(Z/k)  '   =  1.0. 

Figure  18  shows  the  graph  for  the  experiment  conducted  on 

20  April.   The  bottoming  out  occurs  for  the  aperture  size 

1/2 
nearest,  but  not  less  than,  D/(Z/k)  '   =  1.0.   Note,  for 

1/2 
this  path  length  (Z/k)    -5.8  mm,  and  although  the  5  mm 

1/2 
aperture  size  is  nearer  to  (Z/k)  '   than  the  7  mm  size, 

previous  runs  have  shown  that  the  behavior  of  G(D)  below 

1/2 
D/(Z/k)     =  1.0  is  erratic.   Thus,  the  7  mm  size  is  nearer 

to  the  minimum  G(D)  value  that  corresponds  to  the  bottoming 

1/2 
out  point  at  D  =  (Z/k)     than  is  the  5  mm  size.   As  a 

further  comparison,  the  corridor  runs  conducted  on 

7  February  and  (7-9)  April  at  a  path  length  of  280  m  were 

plotted  with  this  run  of  20  April  at  the  path  length  of  336  m, 

shown  in  Figure  19.   It  appears  that  the  bottoming  out  point 

1/2 
occurs  very  close  to  D/(Z/k)     =1.0. 

[G(D)]  -n  for  this  experiment,  with  respect  to 
a  20  mm  D,  was  .666  and  the  average  a.,  was 

^o-rtFCD  ,t)]>  =  .162.   This  does  not  fit  in  with  the  results 
given  in  Table  III;  however,  the  propagation  path  and  turbu- 
lence conditions  were  different.   To  verify  a  predictable 

increase  of  [G(D)1    with  <onT>        would  require  repetition 
L  v  J  Jmin        iol 

of  the  experiment  under  identical  experimental  conditions 
and  over  the  same  propagation  path  where  the  turbulence 
would  have  the  same  characteristics.   Time  considerations 
precluded  this  procedure. 
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d.   Over  Water  (Land  to  Land) 

(1)  2220,  5  March  -  0345,  6  March  1976  (Spanagel 
Roof  to  Coast  Guard  Parking  Lot) 

Due  to  the  irregular  propagation  path  used 
and  corresponding  turbulence  conditions  over  a  long  distance, 
this  entire  sequence  was  conducted  under  saturation  con- 
ditions for  a£I  .   <aAI[F(D  ,t)]>  was  1.01  and  [G(D)]  . 
was  .8808  with  respect  to  a  D  of  24.3  mm.   The  graph  of 

"G(D)  vs.  Aperture  Diameter,  DM  is  shown  in  Figure  20. 

1/2 
For  this  path  length,  (Z/k)  '   =  15.3  mm.   The  main  item 

of  interest  for  this  run  is  that,  under  saturation  conditions, 

the  maximum  difference  in  G(D)  values  observed  was  only 

[AG(D)]    =  .167,  which  included  comparisons  of  aperture 

sizes  of  24.3,  20,  16,  12,  7,  6,  5,  4,  3  and  2  millimeters. 

This  indicates  that  aperture  averaging  has  only  a  small 

effect  on  the  measured  signal  variance  under  saturation 

conditions  for  aperture  sizes  less  than  24.3  mm. 

(2)  1550-1640,  11  March  1976  (Monterey  Holiday 
Inn  to  Hopkins  Marine  Laboratory) 

Due  to  the  limited  time  available  for 

conducting  the  experiment  on  this  day,  only  aperture  sizes 

of  41.7,  24.3  and  14.0  millimeters  were  examined.   A  plot 

of  MG(D)  vs.  Aperture  Diameter,  D"  is  shown  in  Figure  21. 

1/2 
The  path  length  was  4230  m,  giving  a  (Z/k)  '      of  20.6  mm. 

Figure  21  contains  insufficient  data  to  draw  definite 

conclusions.   The  same  data  is  plotted  again  in  Figure  24 

in  terms  of  the  normalized  diameter,  and  will  be  discussed 
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at  that  point.   CGCD)]min  and  <a^i fFCD  »t) ]>  for  this  run 
should  not  be  compared  with  other  experiments  in  this  sec- 
tion since  a  different  source  and  fixed  aperture  detector 
system  was  used  and  a.,  is  wavelength  dependent  for  identical 
turbulence  conditions. 

3.   Propagation  in  the  Marine  Boundary  Layer  (R/V  ACANIA 
Experiments; 

a.   1840-1930,  30  March  1976 

The  winds  during  the  time  of  this  experiment 
varied  from  8  to  12  m/sec  in  a  direction  of  about  10  degrees 
off  the  normal  to  the  propagation  path.   The  sky  was  partly 
cloudy,  temperature  about  7°C,  sea  conditions  were  moderate 
with  average  wave  heights  of  about  1  meter,  and  swell  heights 
of  about  2  to  3  meters  moving  in  a  direction  normal  to  the 
propagation  path.   The  motion  of  the  R/V  ACANIA  at  anchor 
made  maintenance  of  gyro  lock-on  difficult.   This  undoubtedly 
caused  intensity  variances  to  be  greater  than  would  be  ob- 
served normally.   The  path  length  on  this  day  was  1600  m. 
The  detectors  were  located  about  7  meters  and  the  source 
about  5  meters  above  the  ocean  surface.   Aperture  sizes 
compared  were  38,  31,  24.3,  17  and  10  millimeters.   For 
this  path  length,  (Z/k)1^2  =  12.7  mm.   The  MG(D)  vs. 
Aperture  Diameter,  D"  curve  is  shown  in  Figure  22.   This 
is  not  typical  of  past  aperture  averaging  results,  but 
resembles  those  obtained  under  saturated  conditions.   How- 
ever, <a,  -r  [F(D  ,  t)  ]  >  for  this  experiment  was  .149,  well 
below  the  saturation  region.   It  is  felt  that  the  sea 
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conditions,  rapid  ship  movements,  and  inability  of  the  gyro- 
scope lock-on  system  to  work  effectively  in  moderate  seas 
may  have  obscured  any  detectable  aperture  averaging  effects. 

b.  1340-1415,  2  April  1976 

Wind  and  sea  conditions  on  this  day  were 
identical  to  those  of  30  March.   The  sky  being  clear  was 
the  only  difference  in  weather  conditions.   The  sea  state 
was  still  moderate  and  there  was  still  considerable  ship 
motion.   The  path  length  on  this  day  was  1600  m.   The  gyro- 
scope lock-on  system  had  been  improved  slightly  by  readjusting 
the  existing  circuitry;  all  other  experimental  conditions 
were  as  described  for  30  March.   Aperture  sizes  used  for 
•  this  run  were  40,  30,  20  and  10  mm.   Figure  23  shows 

MG(D)  vs.  Aperture  Diameter,  D"  and  Figure  24  compares 

1/2 
MG(D)  vs.  D/(Z/k)  '    "  for  this  run  with  the  data  taken  over 

water  with  a  path  length  of  4229  m  on  11  March.   Once  again, 

1/2 
it  appears  the  normalization  factor  (Z/k)     successfully 

predicts  the  diameter  of  the  aperture  that  best  approximates 

a  vanishingly  small  receiver. 

c.  1445-1657,  27  April  1976 

The  sky  was  clear  at  1445  but  became  overcast 
at  about  1545.   Winds  were  8  to  10  m/sec  blowing  normal 
to  the  propagation  path;  waves  also  moved  normal  to  the  path 
with  a  height  of  about  1  meter  and  swells  about  2  meters 
high.   The  ship  rolled  considerably  in  the  swells  and  once 
again  gyroscope  lock-on  was  difficult.   The  ship  was 
anchored,  path  length  was  1650  m,  and  aperture  diameters 
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examined  were  40,  35,  30,  25,  20,  16,  13,  10,  5  and  4 
millimeters.   "G(D)  vs.  D/(Z/k)  '    "   for  this  run  is  shown  in 
Figure  25.   G(D)  values  seem  erratic;  this  is  attributed  to 
difficulty  in  stabilizing  the  detector  and  tracking  the 
ACANIA  from  shore. 

d.   1830-2024,  29  April  1976 

The  sky  was  clear,  wind  speed  throughout  the 
experiment  was  3  to  4  m/sec,  and  path  length  was  reduced 
to  1442  m.   Sea  conditions  improved  slightly  with  wave 
heights  down  to  about  .7  m  and  swells  of  about  1.5  m.   The 
wave  direction  was  about  15  degrees  off  the  normal  to  the 
propagation  path  as  was  the  wind.   Aperture  sizes  of  40,  35 

30,  25,  20,  16  and  13  millimeters  were  used.   MG(D)  vs. 

1/2 
D/(Z/k)  *  "  is  plotted  in  Figure  26.   The  35  mm  aperture 

size  was  used  as  D   for  the  G(D)  calculations  since  it 

appeared  to  give  more  credible  overall  values  based  on 

past  experiments.   The  G(D)  value  for  the  40  mm  aperture 

is  plotted  relative  to  the  35  mm  aperture  as  D  .   With  the 

exception  of  the  40  mm  aperture  [ j-y^-   =  3.33]  and  the 

(Z/k)1/Z 

16  mm  aperture  [ -r-j-™-   =  1.33],  the  curve  appears  to 

(Z/k)i/Z 

parallel  results  of  past  runs,  i.e.,  reaches  a  minimum  at 

-r-r-yy  ~    1.0  and  shows  a  nearly  linear  decrease  of  G(D) 

(Z/k)i/Z 

with  decreasing  aperture  diameter  to  this  point.   For  this 

run,  [G(D)]min  =  .721  and  <a ^  [F (DQ , t) ] >  =  .145. 
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e.   1044-1311,  30  April  1976 

The  sea  state  was  relatively  calm  with  wave 
heights  of  about  .5  m  and  no  apparent  swells.   The  sky  was 
clear,  temperature  about  21°C  and  winds  were  3  to  4  m/sec. 
The  calm  sea  state  enabled  the  gyroscope  tracking  system  to 
work  well  in  the  automatic  mode  and  results  more  consistent 
with  past  sequences  were  obtained  on  this  day.   The  ship  was 
at  anchor  and  the  propagation  path  length  was  1625  m. 

Aperture  sizes  examined  were  40,  16,  13,  10  and  5  milli- 

1/2 
meters  and  (Z/k)     for  this  path  length  was  12.8  mm. 

The  plot  is  shown  in  Figure  27.   Bottoming  out  occurs  near 

D/(Z/k)1/2  =  1.0,  [G(D)]min  =  .766,  and  for  this  run 

«7AI[F(D 0,t)]>  =  .264.   A  comparison  of  [G(D)]m;Ln  and 

<Oqt>   for  the  three  R/V  ACANIA  experiments  is  given  in 

Table  IV. 

TABLE  IV 

D^e  [G(D)]min  <a£I[F(D0,t)]> 

2  Apr  .792  .726 

30  Apr  .766  .264 

29  Apr  .721  .145 

From  the  data  of  Table  IV  it  appears  again  that  in  general 

[G(D)1™-;   increases  with  <o0T>    . 
L  v  *  Jmin  3ii 
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4.   Comparison  of  Theoretical  Aperture  Averaging  Factors 
with  hxpenmental  Results- 

a.   Lutomirski  and  Yura  Model  [Ref.  1] 

For  aperture  diameters  on  the  order  of  the 

intensity  correlation  length,  as  utilized  during  these 

experiments,  Lutomirski  and  Yura  "recommend"  an  aperture 

averaging  factor  of 


(46)  A  =  1 7  . 


Table  V  lists  Lutomirski  and  Yura's  theoretical  value  of  A 

computed  for  our  experiments  and  the  observed  value  of  A. 

2     2 
A,  theoretically  defined  as  a     /o  ,  is  equivalent  to 

so 

our  [G(D)1  •  /[G(D  )]  where  D^  replaces  D  in  equation  (46). 
L  v  y  Jmm  LVo         o   ^ 

The  theoretical  values  do  not  correlate  well  with  the 

empirical  values.   Further,  Lutomirski  and  Yura's  theoretical 

model  does  not  contain  any  turbulence  dependence  which  appears 

to  be  necessary  in  view  of  the  apparent  dependence  of 

[G(D)1  •   on  <anT>  seen  empirically. 
L  *■  J  Jmin      11  r       / 

b.   Fried's  Model  [Ref.  21] 

To  compare  empirical  results  with  Fried's 
theoretically  derived  plot  of  "0  vs.  D/(4Z/k)  '    " 

(Figure  4) ,  whose  curves  are  based  on  an  aperture  diameter 

1/2 
of  D  =  .l(4Z/k)  '   as  representing  a  vanishingly  small 

receiver,  it  was  necessary  to  use  the  experimental  aperture 
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size  for  each  run  that  gave  the  minimum  G(D)  value  as  the 
vanishingly  small  aperture  diameter.   The  effect  of  this  is 
to  shift  the  empirical  curves  to  the  right  of  the  theoretical 
curve.   Figures  28  and  29  show  Fried' s  theoretical  plot 
for  log  amplitude  variance  C„(0)  =  2~   as  compared  with 
experimental  values  for  various  runs.   For  this  graph, 

aperture  diameters  were  appropriately  normalized  for  com- 

1/2 
parison  by  the  (4Z/k)     factor.   Using  equation  (45), 

0   was  calculated  as 


(47) 


gsW,tQ)] 
°sZfV(D',t0)]  ' 


2 
where  a   [V(D',t  )]  is  the  signal  variance  measured  at 

time  t   and  with  variable  aperture  diameter  D'   which 

o  c  * 

corresponds  to  the  aperture  diameter  that  best  approximated 
a  vanishingly  small  receiver,  i.e.,  corresponded  to 
[G(D)]m;in-   as  [V(D,t  J]  is  the  signal  variance  that  would 

have  been  measured  at  time  t  had  aperture  size  D  been  used, 

o      r 

and  was  found  from  the  relation 

7  9  o    2[F(D  ,t  )] 

(48)     a s  [V(D,tJ]  =  a  Z[V(D,t)]  •  -J- °   °    , 

S  asZ[F(DQ,t)] 


where  again,  F  is  the  fixed  aperture  detector. 

Empirical  results  appear  to  follow  Fried' s 
theoretical  model  generally  well.   In  contrast  with 


68 


Lutomirski  and  Yura's  aperture  averaging  model,  Fried 
acknowledges  the  dependence  of  the  aperture  averaging  factor 
on  the  turbulence  conditions,  as  represented  by  different 
curves  for  particular  values  of  Cg(0)  in  Figure  4.   Fried's 
model  indicates  that  as  C.,(0)  increases,  6  increases  for  a 
specific  aperture  diameter,  or  the  inverse  of  the  aperture 
averaging  factor  decreases.   This  is  in  agreement  with 
empirical  results  noted  in  Sections  2. a. (3)  and  3.e.  above. 
From  Figure  28  it  is  evident  that  the  aperture  averaging 
curve  for  7  February  lies  below  that  of  10  February  and  above 
that  of  (7-9)  April.   The   0  values  are  therefore  in 

general  larger  for  10  February  than  for  7  February,  i.e., 

2     2 
the  ratio  of  a     /a    is  larger,  and  the  inverse  of  the 

So 
aperture  averaging  factor  is  less,  as  predicted  by  Fried's 
theory. 

B.   APERTURE  AVERAGING  EFFECTS  ON  THE  RELATIVE  TEMPORAL - 
FREQUENCY  POWER  SPECTRUM 

1 .   General  Results 

The  spectrum  analyzer  frequency  limits  were  selected 

to  examine  frequencies  within  the  range  10  Hz  to  2.56  KHz. 

The  demodulator  in  the  signal  processing  system  eliminates 

frequencies  above  2  KHz,  coinciding  with  the  original 

estimate  of  the  Fan-Scan  frequency.   The  bandwidth  is  further 

reduced  to  a  maximum  frequency  of  1  KHz  for  experiments 

conducted  using  the  Sanborn  HP- 3914  Recorder.   Power 

spectra  from  all  frequency  analysis  experiments  show  the 
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actual  Fan-Scan  frequency  to  be  1.92  KHz,  clearly  evident 
as  a  noise  spike  in  Figures  34,  35  and  36.   Previous  work 
at  NPS  by  Schroeder  [Ref.  28]  examined  frequencies  from 
20  Hz  to  1  KHz. 

Results  of  the  experiments  conducted  on  (7-9)  April, 
20  April  and  (29-30)  April,  under  varying  turbulence 
conditions,  all  indicate  that  aperture  averaging  has  a 
discernable  effect  on  the  observed  power  spectrum.   Empiri- 
cal results  support  the  theoretical  prediction  that  a 
corresponding  reduction  in  high  frequency  power  results  as 
aperture  diameter  is  increased.   The  reduction  in  power, 
however,  appears  to  occur  over  a  large  portion  of  the  fre- 
quency range,  not  solely  on  the  high  frequency  end.   Figures 
30,  31,  32  and  33  show  the  effect  of  aperture  averaging  on 
the  relative  power  spectrum. 

2.   Turbulence-Dependent  Forms  of  the  Frequency  Spectrum 
As  previously  remarked,  the  three  frequency  analysis 
experiments  were  conducted  at  different  locations  and  under 
different  turbulence  conditions.   The  resulting  power 
spectra  for  each  run  were  similarly  different.   Typical 
spectra  for  the  corridor  experiment  (inhomogeneous  turbulence) , 
the  Navy  Annex  run  (over  land  in  homogeneous  turbulence) ,  and 
the  R/V  ACANIA  data  (over  water  in  homogeneous  turbulence) , 
are  shown  in  Figures  34,  35  and  36,  respectively.   The 
significant  item  of  note  in  Figure  34  is  the  uniformity  of 
the  spectrum  below  .85  KHz.   The  prevalence  of  noise  spikes 
at  110  Hz  and  120  Hz  intervals  in  Figure  35  is  believed  due 
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to  the  extremely  noisy,  gas-powered  portable  generator 
used  to  provide  power  to  the  electronic  equipment  during 
this  experiment.   Also  peculiar  to  this  experiment  were 
characteristically  low  magnitude  power  spectra.   Figure  36 
shows  an  approximately  linear  decrease  in  power  with 
increasing  frequency. 

Eliminating  the  noise  spikes  in  Figure  35,  it 
appears  that  the  spectrum  assumes  the  same  general  form  as 
that  of  Figure  34.   A  significant  difference  in  the  experi- 
mental conditions  for  9  April  and  20  April  as  opposed  to 
those  of  30  April  was  the  average  wind  speed  present. 
On  9  April  and  20  April  winds  were  negligible.   On  30  April, 
conducted  at  sea  on  board  the  ACANIA,  average  wind  speeds 
were  about  4  m/sec.   Further  experiments  need  to  be  con- 
ducted to  examine  the  effect  of  wind  speed  on  the  shape  of 
the  power  spectrum. 

3.   Comparison  with  Clifford's  Theoretical  Model 
Clifford's  theoretical  model  for  the  temporal- 
frequency  power  spectrum  for  spherical  and  plane  waves 
are  given  by  equations  (32)  through  (35)  .   The  values  of  f 
for  9  April,  20  April  and  30  April,  based  upon  respective 
average  wind  speeds  of  1  m/sec,  1.5  m/sec  and  4  m/sec  are 
30  Hz,  41  Hz  and  53  Hz,  respectively.   For  comparison  with 
the  theoretical  model  (Figure  3) ,  typical  spectra  for  these 
experiments  were  plotted  vs.  "Log  ft"  and  are  shown  in 
Figure  37.   The  vertical  axis  is  the  percentage  of  power 
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at  a  particular  frequency  compared  to  the  maximum  power 
attainable  as  determined  by  the  spectrum  analyzer's 
"Y-axis"  voltage  setting. 

There  appear  to  be  two  significant  differences 
between  the  empirical  results  of  Figure  37  and  Clifford's 
theoretical  spectra  of  Figure  3.   First,  Clifford  indicates 
that  the  rapid  downward  slope  of  the  spectrum  begins  at 
log  0,   ~-0.1  for  plane  waves  and  log  0,   ~  0 . 2  for  spherical 
waves.   Empirical  results  show  the  drop  occurring  at 
log  ft  ~  3.0  for  the  (7-9)  April  and  20  April  experiments. 
The  spectrum  of  (29-30)  April  shows  an  overall  negative 

slope  for  all  frequencies.   Secondly,  the  slope  appears  to 

-8/3 
decrease  less  rapidly  than  the  theoretical  (f)   '  .   The 

slope  for  (7-9)  April  is  of  the  order  of  (f)"  '    and  for 

20  April  and  (29-30)  April,  (f)~2/13.   Further  data 

collection  is  necessary  to  examine  if,  and  to  what  extent, 

various  turbulence  parameters  affect  the  magnitude  of  the 

negative  slope.   It  does  appear,  however,  that  Clifford's 

-8/3 
(f)     is  too  rapid  a  decrease. 

4.   Differences  in  Magnitude  of  the  Aperture  Averaging 
Ettect 

Examination  of  Figures  30  through  33  indicates 

that  there  are  significant  differences  between  aperture 

averaging  effects  on  the  relative  power  spectrum  as  seen 

under  the  various  experimental  conditions.   In  an  attempt 

to  determine  a  factor  common  to  the  variances  and  to  the 

power  reductions  which  were  not  directly  proportional  to 
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aperture  diameter,  parameters  such  as  wind  speed,  <^j?T>  , 
CN,  a. y  and  G(D)  were  examined.   No  factor  was  found 
that  individually  accounts  for  the  variations.   A  detailed 
analysis  of  all  factors  concurrently  was  not  undertaken 
due  to  time  limitations. 
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V.   SUMMARY  AND  CONCLUSIONS 

A.   APERTURE  AVERAGING  EFFECTS  ON  SCINTILLATION 

The  variance  of  log -amplitude  fluctuations  that  would 
theoretically  be  observed  if  using  a  vanishingly  small 

receiver  can  best  be  approximated  by  an  aperture  of  diameter 

1/2 

(Z/k)    .   A  diameter  of  this  size  eliminates  aperture 

averaging  effects  present  when  using  apertures  larger  than 

1/2 
(Z/k)    ,  and  that  may  be  present  when  using  aperture 

1/2 
diameters  smaller  than  (Z/k)    .   Experiments  conducted 

in  an  enclosed  laboratory  in  inhomogeneous  turbulence  and 

field  experiments  conducted  in  homogeneous  turbulence, 

both  over  land  and  in  the  marine  boundary  layer,  support 

this  prediction  from  the  theoretical  model  of  Lutomirski 

and  Yura  [Ref.  1] . 

Fried's  theoretical  model  of  aperture  averaging 

[Ref.  21]  predicts  that  0,  the  aperture  averaging  factor, 

increases  continuously  as  aperture  size  is  decreased  from 

an  arbitrarily  large  aperture  diameter.   Empirical  results 

indicate  that  the  aperture  averaging  factor  does  not  increase 

continuously  as  aperture  size  is  reduced  to  a  vanishingly 

small  diameter,  but  increases  continuously  until  an  aperture 

1/2 
diameter  (Z/k)     is  reached.   A  further  reduction  in 

1/2 
aperture  diameter  [below  (Z/k)  '  ]  causes  0  to  (1)  reverse 

direction  and  decrease  or  (2)  remain  near  the  same  value 

of  0. 
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In  general,  the  value  of  the  aperture  averaging  factor 

increases  as  <oJ,T>  increases.   Put  another  way,  for  an 

1/2 
aperture  of  diameter  larger  than  (Z/k)    ,  the  aperture 

averaging  correction  factor  decreases  as  the  amount  of 

turbulence  increases.   Under  weak  turbulence  conditions, 

the  effect  of  aperture  averaging  is  quite  significant.   For 

strong  scintillation  conditions,  the  aperture  averaging 

effect  is  decreased  and  becomes  less  significant. 

The  ratio  G(Df )/ [G(D) ]  .   gives  the  factor  by  which  the 

measured  signal  variance  (or  other  turbulence  parameter) , 

1/2 
measured  with  an  aperture  of  diameter  D',  D'  >  (Z/k)  '  , 

must  be  multiplied  to  represent  the  value  that  would  have 

been  measured  had  a  vanishingly  small  receiver  been  used. 

It  is  important  to  note  that  the  value  of  [G(D)]  .   is  a 

function  of  the  largest  aperture  diameter  used  (D  ) ,  and 

should  not  be  viewed  as  an  absolute  value  for  a  particular 

set  of  turbulence  parameters.   In  general,  the  value  of 

fG(D)]  .   decreases  as  D   is  increased. 
L    Jmin  o 

B.   APERTURE  AVERAGING  EFFECTS  ON  THE  TEMPORAL -FREQUENCY 
POWER  SPECTRUM 

An  increase  in  aperture  diameter  appears  to  decrease 

the  magnitude  of  power  in  the  frequency  range 

.5  KHz  <  f  <  2  KHz  of  the  temporal -frequency  power 

spectrum  for  scintillation  measurements.   The  difference 

in  power  at  a  particular  frequency  for  two  given  aperture 

diameters  differed  significantly  for  different  experimental 

and  turbulence  conditions.   The  reason  for  the  disparity 

was  not  determined  due  to  time  considerations. 
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Spectra  obtained  from  experiments  in  the  corridor 
laboratory,  over  land,  and  in  the  marine  boundary  surface 
layer  indicate  a  negative  sloped  decrease  in  power  at 

higher  frequencies  for  the  temporal- frequency  power 

-  8/3 
spectrum  less  rapid  than  the  (f)   '   predicted  by 

Clifford  [Ref.  19].   Further  investigation  is  required 

to  determine  the  effect  changes  in  certain  turbulence 

parameters  have  on  the  form  of  the  power  spectrum. 

C.   EXTENSION  OF  PROJECT 

Comparisons  of  [G(D)1  .   with  <onT>   indicated  that, 
r  L     Jmin        Jc  1 

for  a  given  aperture  size  larger  than  that  which  approxi- 
mates a  vanishingly  small  receiver,  the  aperture  averaging 
correction  factor  is  not  a  constant,  but  is  dependent  upon 
the  level  of  turbulence  present.   Further  experimentation 

is  needed  to  compare  [G(D)1    values,  for  a  given  maximum 

c  L     Jmm 

reference  diameter  D  ,  with  <^5t>  values,  to  determine  the 
relationship  between  the  level  of  turbulence  and  the  aperture 
averaging  correction  factor  for  a  given  aperture  diameter  D  . 

It  seems  clear  from  this  report  and  other  published 
data  that  aperture  averaging  reduces  the  high  frequency 
power  content  in  the  temporal -frequency  power  spectrum 
and  eliminates  certain  high  frequencies  proportional  to 
increasing  aperture  diameter.   However,  experiments  conducted 
under  different  turbulence  conditions  and  over  different 
propagation  paths  indicate  that  the  power  spectrum  assumes 
different  general  shapes.   This  needs  further  investigation 
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through  the  examination  of  power  spectra  obtained  over  the 
same  propagation  path  under  various  turbulence  conditions, 
and  under  the  same  turbulence  conditions  over  various 
propagation  paths.   In  conjunction  with  this,  Clifford's 
theoretical  model  of  the  temporal -frequency  power  spectrum 
for  plane  and  spherical  waves  should  be  examined  for  validity 
more  thoroughly  through  empirical  study. 
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Figure  3 
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Clifford's  theoretical  temporal  power  spectrum  of  log 
amplitude  fluctuations  for  plane-  and  spherical-wave 
propagation. 
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Figure  4 
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Figure  8 
Electronic  Components:   Aperture  Averaging  Experiment 
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APPENDIX  B 

TEST  PROCEDURE  AND  PERFORMANCE  CHARACTERISTICS 
OF  THE  TEXAS  INSTRUMENTS  TYPE  TIXL69  LARGE-AREA 
SILICON  AVALANCHE  PHOTODIODE 


A  series  of  tests  was  conducted  to  evaluate  the 
detection  capabilities  and  performance  characteristics  of 
the  TIXL69  avalanche  photodiode  detector  and  to  compare 
them  with  that  of  the  PIN-5DP  photodiode  detector.   The 
two  detectors  were  evaluated  at  a  signal  intensity  compara- 
ble to  that  anticipated  during  a  typical  at-sea  data  run, 

-9     2 
i.e.,  about  1.4X10   W/cm   (the  intensity  of  a  5  mW  He-Ne 

o 

(6328A)  laser  source  at  a  distance  of  about  .8  nautical 
miles  and  with  a  beam  divergence  of  7  milliradians) . 

The  circuit  used  with  the  avalanche  detector  is 
illustrated  in  Figure  Bl(a).   The  reverse  bias  voltage 
source  was  a  Hewlett  Packard  Harrison  6515A  DC  Power 
Supply. 

Dark  current  as  a  function  of  reverse  bias  voltage 

was  measured  first  to  see  if  it  compared  with  the  manu- 

-9  -7 

facturer's  stated  values  of  3.5X10  A  minimum  to  1.1X10  A 

maximum.   The  test  circuit  used  is  shown  in  Figure  Bl(b). 

During  this  experiment  the  avalanche  detector  was  covered 

by  black  cloth  to  keep  light  from  hitting  the  detector 

surface.   There  were  no  windows  in  the  enclosed  room,  all 

room  lights  were  out,  and  the  room  was  completely  dark. 

A  flashlight  was  used  for  recording  data.   To  ensure  that 
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the  voltage  applied  was  actually  occurring  as  a  voltage 
drop  across  the  avalanche  photodiode,  a  Hewlett  Packard 
419A  Null  Voltmeter  was  connected  to  the  circuit  as  shown. 
The  voltmeter  indicated  that  the  avalanche  photodiode  was 
indeed  receiving  the  full  potential  drop.   Using  a  Keithley 
Instruments  600B  Electrometer  to  measure  the  current,  it 
appeared  that  values  in  the  nanoamp  region  were  occurring. 
However,  the  current  indicator  was  not  stable  and  repetitive 
current  readings  at  the  same  reverse  bias  voltage  fluctuated 
greatly.   It  was  decided  that  a  more  accurate  ammeter  was 
needed  and  a  Keithley  Instruments  410  Micro-Microammeter  was 
obtained.   Steady,  repeatable  values  of  current  were  now 
attained  for  specific  reverse  bias  voltages. 

As  shown  in  Figure  B2,  dark  current  was  linear  with 

applied  reverse  bias  voltage  from  20  volts  through  157.5 

3 
volts,  then  increased  by  a  factor  of  about  10   at  an 

applied  voltage  of  160  volts.   Dark  current  values  were 
higher  than  that  stated  by  the  manufacturer  by  a  factor  of 
10   at  low  applied  reverse  bias  voltages,  decreasing  to  a 
factor  of  10  higher  at  high  values.   For  this  particular 
detector,  breakdown  occurs  at  a  reverse  bias  voltage  of 
160  volts  vice  the  manufacturer's  indicated  166  volts. 

The  HP  Harrison  6515A  DC  voltage  source  was  then  re- 
placed by  a  battery  box  consisting  of  a  series  of  seven 
22.5  volt  and  one  1.5  volt  batteries  to  be  used  with  the 
avalanche  detector  during  actual  aperture  averaging 
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experiments.   The  "Dark  Current  vs.  Reverse  Bias  Voltage" 
experiment  was  repeated.   The  same  equipment  and  circuit 
as  described  above  was  used  and  dark  current  values  obtained 
with  the  battery  box  were  identical  to  those  obtained  with 
the  HP  Harrison  6515A  DC  source. 

With  the  breakdown  voltage  known,  the  next  test  was  to 
examine  the  "Signal-to-Noise  Ratio  vs.  Reverse  Bias  Voltage" 
at  signal  intensities  similar  to  those  previously  described 
under  at-sea  conditions.   The  experimental  setup  was  as 
shown  in  Figure  B3.   The  signal  source  was  a  . 5  mW  He-Ne 

o 

(6328A)  laser.   The  laser  beam  was  passed  through  a  +8  cm 
focal  length  lens  followed  by  a  signal  chopper  located  just 
beyond  the  focal  point  of  the  lens,  then  proceeded  a  distance 
of  3.35  meters  to  the  avalanche  detector.   The  avalanche 
detector  was  mounted  on  a  stand  that  permitted  horizontal  and 
vertical  movement  in  order  to  maximize  signal  intensity 
appearing  on  the  oscilloscope.   The  signal  from  the  detector 
was  sent  to  a  PAR  (Princeton  Applied  Research)  Model  113 
Pre-Amp,  then  to  a  Monsanto  AM-6419/USM-368  oscilloscope 
for  display.   External  triggering  was  used  on  the  oscillo- 
scope.  To  obtain  the  desired  signal  intensity,  a  Jodon 
Model  PM-550  Digital  Optical  Power  Monitor,  placed  at  the 
location  of  the  avalanche  detector,  measured  signal  intensity 
directly.   With  the  uninterrupted  laser  beam  incident  on  the 
Jodon  detector  at  close  range  (about  10  cm) ,  the  Power 
Monitor  indicated  a  maximum  power  output  of  .39  mW.   The 
diameter  of  the  laser  beam  was  6  mm  wide  while  the  diameter 
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of  the  active  detection  area  of  the  Jodon  detector  is  3  mm. 
Since  the  laser  beam  is  theoretically  (and  appeared  to  be 
from  visual  inspection)  Gaussian,  roughly  87%  of  the  output 
intensity  was  incident  on  the  Jodon  detector  active  area, 

covering  it  completely,  and  giving  a  source  power  output 

-3    2 
of  5.5X10   W/cm  .   At  the  position  of  the  avalanche  detector, 

the  intensity  had  decreased,  due  to  beam  expansion  from  the 

-  3     2 
unfocused  optics  of  the  source  output,  to  1.0X10  W/cm  . 

With  the  chopper  added  and  running  at  normal  operating 

-4    2 
speed,  the  intensity  reduced  further  to  5.7X10  W/cm  . 

With  the  diverging  lens  in  place,  to  complete  the  setup, 

_  5 
intensity  decreased  to  3X10   W/cm  .   To  reduce  the  intensity 

further,  neutral  density  filters  (Kodak  Wrattan  Filter  (in 

B  glass)  N.D.4,  N.D.3  and  N.D.2)  were  used  singly  and  in 

combination.    With  the  N.D.4  filter  inserted  (reduces 

-4 
intensity  by  10   )  between  the  laser  and  the  +8  cm  lens, 

-  9     2 
signal  intensity  was  reduced  to  3X10  W/cm  ,  as  desired. 

When  using  the  PAR,  it  was  found  that  a  saturated 
signal  occurred  on  the  oscilloscope  if  PAR  Gain  setting 
was  too  high.   A  saturated  signal  was  identifiable  as  a 
rectangular-shaped  8.75  volt  signal  with  normally  visible 
scintillation  absent. 

On  26  January  1976  an  experiment  comparing  "S/N  vs. 
Reverse  Bias  Voltage"  at  a  given  signal  intensity  for  the 
avalanche  detector  was  conducted  followed  by  a  comparison 
of  "S/N  of  Avalanche  Detector  vs.  S/N  of  PIN-5DP  Detector" 
at  varying  intensities.   In  both  experiments  the  equipment 
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setup  of  Figure  B3  was  used.   PAR  settings  were:   L-F 

ROLL  OFF--100,  H-F  ROLL  OFF--300K,  Gain—Variable.   The 

Harrison  6515A  DC  Power  Supply  was  used  to  supply  reverse 

bias  voltage  to  the  avalanche  detector  and  a  6328A  filter 

was  used  in  front  of  the  avalanche  and  PIN  photodiode 

detectors.   No  lens  system  was  used  to  focus  the  incoming 

light  onto  the  detector  active  area  for  either  detector. 

Table  BI  gives  data  collected  during  the  "S/N  vs.  Reverse 

Bias  Voltage"  for  the  avalanche  detector,  shown  graphically 

in  Figure  B4.   In  Table  BI ,  note  (1)  the  lower  PAR  Gain 

settings  required  to  avoid  saturation  as  V,  .    increases 

d  i  as 

beyond  120  volts,  and  (2)  the  increase  in  noise  voltage 

as  VK •    increases.   From  Figure  B4  it  can  be  seen  that  S/N 
d  x  as 

increases  fairly  linearly  from  a  reverse  bias  voltage  of 
80  volts  through  159  volts.   Table  BII  is  a  comparison  of 
the  performance  of  the  PIN  photodiode  detector  and  the 
avalanche  photodiode  detector.   Both  detectors  were  located 
at  the  same  distance  from  the  signal  source  and  both  used 
the  same  PAR  and  oscilloscope  for  signal  output.   It  is 
evident  that  the  avalanche  detector,  when  set  at  its 
optimum  reverse  bias  voltage,  gives  a  significant  improve- 
ment in  performance  over  the  PIN  photodiode  detector, 
especially  at  the  weaker  intensities  characteristic  of 
actual  field  experiments. 

A  responsivity  curve  was  plotted  to  compare  experimental 
results  with  the  manufacturer's  stated  performance  charac- 
teristics.  From  Figure  B5  it  can  be  seen  that  experimental 
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results  follow  very  closely  the  predicted  values  with  the 
following  exceptions:   (1)  actual  signal  gains  corresponding 
to  manufacturer's  highest  predictions  were  not  attained  and 
(2)  the  maximum  reverse  bias  voltage  that  could  be  applied 
was  less  than  manufacturer  indications. 

In  summary,  the  TIXL69  avalanche  detector  proved  to 
give  a  significant  improvement  over  the  PIN-5DP  photo- 
diode  detector  for  low  intensity  signal  detection,  as 
anticipated.   The  maximum  (and  optimum),  reverse  bias  voltage 
that  can  be  applied  is  159  volts,  slightly  less  than  previous 
ly  expected.   Measurements  of  responsivity  agree  with 
manufacturer's  stated  values  but  do  not  attain  the  highest 
values  predicted. 
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Figure  Bl(a) 
Circuit  Diagram:   Avalanche  Photodiode  Detector 
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Figure  Bl(b) 
Test  Setup:   Dark  Current  Experiment 
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TABLE  BI 
S/N  VS.  VBIAS 
'FOR  AVALANCHE  DETECTOR 


26  January  1976,  Source:   .5mW,  6328  A  He-Ne  laser. 
Using  chopped  signal  and  diverging  lens,  as  shown  in 
Figure  B3,  signal  intensity  using  N.D.F  4  ~  3X10_9W/cm2 


VBIAS 

PAR 
GAIN 

V(NOISE) 

V(SIGNAL) 

S/N 

20V 

10K 

.55V 

1.4V 

2.5 

50V 

10K 

.60V 

1.55V 

2.6 

80V 

10K 

.65V 

1.9V 

2.9 

100V 

10K 

.70V 

2.75V 

3.9 

120V 

10K 

1.05V 

5.25V 

5.0 

140V 

5K 

1.25V 

7.0V 

5.6 

150V 

2K 

1.10V 

6.5V 

5.9 

155V 

IK 

1.00V 

6.0V 

6.0 

159V 

200 

.70V 

4.35V 

6.2 

160V 

NO 

SIGNAL 
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Figure  B5 

Comparison  of  Manufacturer's  Predicted  and 
Empirical  Responsivities 


Solid  Line  =  Manufacturer's     ©=  Empirical 
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APPENDIX  C 
EQUIPMENT  PERFORMANCE  CHARACTERISTICS 

1.  Precision  Instruments,  Model  No.  PI-6204,  four-channel 
recorder: 

Input  sensitivity,  .05  to  10  volts  rms ; 

Frequency  response,  DC  to  10  KHz  +3  dB  at  37.5  ips 

on  FM  Record; 
Signal  to  Noise  Ratio,  35  dB  at  37.5  ips,  300  Hz 

to  100  KHz. 

2.  Sanborn  HP-3914,  14-channel  Analog  Tape  Recorder: 

Maximum  Input,  5  volts  (peak  to  peak); 

Gain,  unity  from  0  to  850  Hz,  -.23  dB  at  1  KHz; 

Tape  Speed  7.5  ips. 
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